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Abstract. A multi-line millimeter- wave study of the nearby starburst galaxy NGC 4945 has been carried out using 
the Swedish-ESO Submillimeter Telescope (SEST). The study covers the frequency range from 82 GHz to 354 GHz 
and includes 80 transitions of 19 molecules. 1.3 mm continuum data of the nuclear source are also presented. An 
analysis of CO and 1.3 mm continuum fluxes indicates that the conversion factor between H2 column density 
and CO J=l-0 integrated intensity is smaller than in the galactic disk by factors of 5-10. A large number of 
molecular species indicate the presence of a prominent high density interstellar gas component characterized 
by TiHa ^lO^cm""^. Some spectra show Gaussian profiles. Others exhibit two main velocity components, one at 
~450kms~'^, the other at ~710kms~^. While the gas in the former component has a higher linewidth, the latter 
component arises from gas that is more highly excited as is indicated by HCN, HCO"'' and CN spectra. Abundances 
of molecular species are calculated and compared with abundances observed toward the starburst galaxies NGC 253 
and M 82 and galactic sources. Apparent is an 'overabundance' of HNC in the nuclear environment of NGC 4945. 
While the HNC/HCN J=l-0 line intensity ratio is ~0.5, the HNC/HCN abundance ratio is ~1. From a comparison 
of Ka=0 and 1 HNCO line intensities, an upper limit to the background radiation of 30 K is derived. While HCN 
is subthermally excited (Tox'^SK), CN is even less excited (Tox~3-4K), indicating that it arises from a less dense 
gas component and that its N=2-l line can be optically thin even though its A''=l-0 emission is moderately 
optically thick. Overall, fractional abundances of NGC 4945 suggest that the starburst has reached a stage of 
evolution that is intermediate between those observed in NGC 253 and M82. Carbon, nitrogen, oxygen and sulfur 
isotope ratios are also determined. Within the limits of uncertainty, carbon and oxygen isotope ratios appear to 
be the same in the nuclear regions of NGC 4945 and NGC 253. High ^^0/"0, low "0/^**0 and ^-^N/^^N and 
perhaps also low ^^S/^^S ratios (6.4±0.3, 195±45, 105±25 and 13.5±2.5 in NGC 4945, respectively) appear to be 
characteristic properties of a starburst environment in which massive stars have had sufficient time to affect the 
isotopic composition of the surrounding interstellar medium. 
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1. Introduction 

NGC 4945, a nearby edge-on (i=78°) spiral starburst 
galaxy at a systemic velocity of 14ys ^ 560kms~^ (e.g. 
Ott et al. 2001), is one of the three brightest IRAS 
point sources (^loo^m — 686 Jy; IRAS 1989) beyond the 
Magellanic clouds and a goldmine for molecular cloud re- 
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search. NGC 4945 contains a highly obscured 'compton 
thick' (e.g. Maiolino et al. 1999) Seyfert 2 nucleus and is 
classified as SB(s)cd or SAB(s)cd. Being a member of the 
Centaurus group, its distance was estimated to lie between 
3.9 Mpc (de Vaucouleurs 1964; Bergman et al. 1992) and 
8.1 Mpc (Baan 1985). The value of 6.7Mpc is most often 
used (e.g. Dahlem et al. 1993; Henkel et al. 1994a, here- 
after H94; Ott et al. 2001), but 4Mpc (20" correspond 
to 390 pc) is a more realistic estimation (Bergman et al. 
1992) and will be used throughout the article. 
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Table 1. Observational parameters 



V 






sys 




(GHz) 


(arcsec) 


(pc) 


(K) 




82-99 


63-53 


1220-1030 


150-200 


0.78-0.73 


100-153 


52-34 


1000-660 


300-400 


0.73-0.61 


154-230 


33-22 


640-420 


600-1000 


0.61-0.46 


241-267 


21-19 


410-370 


1500-2300 


0.44-0.41 


330-354 


19-15 


370-290 


-^3000 


0.32-0.30 



Full width to half power (FWHP) beam widths. To estab- 
lish linear seales, L>— 4Mpc was adopted. 

' System temperatures in units of main beam brightness tem- 
perature (Tmb) 

Beam efficiencies derived from measurements of Jupiter (L. 
Knee, priv. comm.) 



Numerous studies from the cm-waveband to hard X- 
rays have been made towards this prominent southern 
galaxy. Associated with the nuclear region is a strong ra- 
dio continuum source with S'i.4ghz=4.8 Jy (Whiteoak & 
Bunton, 1985; Ott et al. 2001) and the first known mega- 
maser (Dos Santos & Lepine 1979; Batchclor et al. 1982), 
detected in the 22 GHz water vapor line. Numerous molec- 
ular emission lines were detected, initially at cm and later 
at mm wavelengths (Whiteoak & Gardner 1975; Whiteoak 
& Wilson 1990; Henkel et al. 1990; H94; Curran et al. 
2001, hereafter COl). An edge-on nuclear gas ring with a 
size of ~ 30" was discovered by Whiteoak ct al. (1990) 
in the J=l-0 line. A detailed analysis of ^'^Q}^'Q> 

and i^ci^O profiles by Bergman et al. (1992) and ad- 
ditional CO observations by Dahlem et al. (1993) and 
Mauersberger et al. (1996a) further support the presence 
of this ring. An HON map indicates that low and high 
density molecular gas coexists in the central 30" (H94). 
On a much smaller scale, the distribution of H2O mega- 
maser emission suggests the presence of a circumnuclear 
disk with a radius of 20 milliarcseconds and a binding 
mass of ~ 10*5 Mq (GreenhiU et al. 1997). 

A nuclear source with a column density of A/'(H) ~ 
4xl0^^cm~^ and possibly disk-like morphology is indi- 
cated by X-ray data (Done et al. 2003). Strongly pro- 
cessed ices were detected near 4.65/im in the more ex- 
tended nuclear environment (Spoon et al. 2003). A conical 
wind-blown cavity was revealed by ground based spectra 
and images at optical and IR wavelengths (Nakai 1989; 
Moorwood ct al. 1996). With the NICMOS Camera on- 
board of the Hubble Space Telescope (HST) a starburst 
ring of size 5" was found in Pa a; the walls of a conical 
cavity blown by supernova driven winds were detected in 
H2 (Marconi et al. 2000). X-ray observations reveal a com- 
plex morphology including a strong variable nuclear source 
at 2-10 kcV and a conically shaped 'plume' at 0.3-2 keV 
(Schurch et al. 2002). 

To better understand the nuclear environment of 
NGC 4945, we present observations of a large number of 
molecular lines towards its central region and calculate 
relative abundances, isotope ratios and H2 densities. 



2. Observations 

2.1. Spectroscopic observations 

The observations were carried out in May 1994, September 
1995, March and .June 1996, January and July 1997, 
January and July 1998, September 2001 and February 
2003 with the 15-m SEST telescope at La Silla. 3 mm 
and 2 mm receivers as well as 3 mm and 1.3 mm receivers 
were employed simultaneously. For the observations at 
330 354 GHz, another SIS receiver was used that was de- 
scribed in detail by Mauersberger et al. (1996a). Two 
acousto-optical backends (AOS) each covered a bandwidth 
of 1 GHz. With a total of 1600 and 1440 channels, the 
channel spacing was 0.68 MHz. 

All measurements were made in a dual beam switch- 
ing mode (switching frequency 6 Hz) with a beam throw 
of 11'40" in azimuth. Since rapid beam switching was 
used in conjunction with reference positions on both sides 
of the galaxy, baselines are of good quality. Calibration 
was obtained with the chopper wheel method (for inter- 
nal consistency and comparisons with other telescopes, 
see e.g. Sect. 3.2 and H94 (their Sect. 2) and Chin et al. 
(1996, their Sect. 2)). Observed frequency ranges, beam 
widths, system temperatures and beam efficiencies are 
given in Table 1. The latter, measured at 94, 115, 230 and 
345 GHz, were interpolated and, if necessary, also extrap- 
olated to convert antenna temperatures (T^) to a main 
beam brightness temperature (T^y,) scale. The pointing 
accuracy, obtained from measurements of the nearby con- 
tinuum source CenA, was mostly better than 10" (see 
Sect. 3.2). 

2.2. 1.3 mm continuum measurements 

In July 1997, we mapped the 1.3 mm continuum emission 
toward the central region of NGC 4945 with a single pixel 
bolometer at the SEST. We first searched for extended 
emission by mapping a 200" x 200" area using the on the 
fly dual beam technique. In addition, the inner 60" x 60" 
were mapped in raster mode yielding a better signal-to- 
noise ratio. 

Dual beam switching was performed using the SEST 
focal plane chopping mirror with an offset of 67.5" in az- 
imuth and 3" in elevation. Pointing was frequently checked 
on Cen A and on NGC 4945 itself. The mean pointing ac- 
curacy was ±2.5". Uranus was used as the primary calibra- 
tor with an assumed flux density of 37Jy/beam. Fitting 
two dimensional (2D) Gaussians to maps of Uranus yields 
a deconvolved full width to half power (FWHP) beam size 
of 27.7". The sensitivity of the bolometer system was de- 
termined to be 190mJy/beams~^. 

The double beam maps were taken in a horizontal sys- 
tem with an extent of 4' x 3', a scanning velocity of 8"/s 
in azimuth, and a step size in elevation of 8". We took 
36 individual maps, which were corrected for atmospheric 
extinction using opacities obtained from skydips, and for 
the SEST gain-elevation curve. Using the MOPSI pack- 
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Fig. 1. CO, CS and N2H"'" spectra measured toward the center of NGC 4945. The velocity resolution for each spectrum 
is ^ 20kms~^, the velocity scale is Local Standard of Rest. Gaussian fits were made towards detected and tentatively 
detected lines, and resulting parameters are listed in Table 2 (same in Figs. 2-5). 



age, the double beam maps were converted into single- 
beam maps and transformed into equatorial coordinates. 
The final coadded map has an rms of 30mJy/boam. The 
raster map was taken in asymmetric on-off mode (off- 
on-on-on, 10 s each phase) with an accumulated integra- 
tion time of 300 s per point, thus reaching an rms of 
llmJy/beam. The 1.3mm emission is slightly more ex- 
tended than the beam size. Fitting a 2D Gaussian to the 
on the fly map the FWHP extent is 27.4" x 30.8" with a 
position angle of 15°, corresponding to a deconvolved size 
of 14" X 20" or 270 pc x 390 pc. Both the on-the-fly map 
and the raster map yield a maximum flux of 1.4 Jy/bcam 
including the line emission. The flux density of the total 
emission mapped in our OTF map is 2.8 ± 0.3 Jy of which 
2.5 ± 0.3 Jy come; from the Gaussian shaped central emis- 
sion. Again these values also include spectral line emission. 
The separation of line and continuum flux and a column 
density estimate are presented in Sect. 8.1. 



3. Results 

3.1. The observed line profiles 

Figs. 1-5 display the spectra from a total of 19 molecu- 
lar species that were measured toward the nuclear region 
of NGC 4945, nominally within 2" of a = 13''05"27.2% 
S = -49°28'05" (J2000.0). According to the overall rota- 
tion curve (see Fig. 5 in Mauersberger et al. 1996a), blue- 
shifted line emission should arise southwest, red-shifted 
emission should originate northeast of the dynamical cen- 
ter. For lineshapes, see H94 (their Sect. 5.2). Gaussian line 
parameters to provide integrated intensities arc given in 
Table 2, including 49 detected line features, 9 tentatively 
detected transitions and 22 undetected lines (for the notes 
to Table 2, see footnote ^) 



^ a) '-|-': detection; ' — ': non-detection; '?': tentative detec- 
tion 'B': blended with other lines, 'N': no baseline. 

b) Integrated from 300 to 800 km s^^ after subtracting a first 
order baseline or a constant offset in T^. The errors were de- 
rived from Gaussian fits and do not account for calibration and 
pointing uncertainties. 

c) Obtained from single component Gaussian fits, vlsr = whel 
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Table 2. Line parameters 
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3.2. Line intensities 

In order to compare our observational data with previ- 
ously published results, we list in Table 3 integrated line 
intensities (/Tmbdw) and peak line temperatures (Tmb)- 
In columns 2 and 5, we display integrated intensities and 
peak line temperatures derived from our data. In columns 
3, 4, 6 and 7, data from COl, H94 and Mauersberger et al. 
(1996a) are given. There is good agreement for the strong 
CO fines, i.e. CO and ^^CO. Differences are 20% or less for 
the J=l-0 and 2-1 transitions. Not surprisingly, for the 
weaker lines, differences are larger. This is likely a conse- 
quence of lower signal-to-noise ratios and a lower pointing 
quality, caused by pointing shifts during one hour long 
integrations and the fact that checks of the shape of the 
profile, showing whether the pointing was stable or not, 
were not possible for individual scans. Because the CO 
and ^"^CO spectra encompass the entire frequency range 
of our observations, we conclude that serious systematic 



— 4.5kms"\ 

d) rms values for a ~20kms~^ channel width on a Tmb scale. 

e) Channel spacings after smoothing as shown in Figs. 1—5. 



calibration problems are absent in H94, Mauersberger et 
al. (1996a), COl and our new data. Significant differences 
in line intensity and lincshapc that arc seen in the weaker 
lines are interpreted as the simultaneous consequence of 
noise and pointing errors and will be analysed in the fol- 
lowing sections. 

3.3. Extent of t lie survey 

The number of measured lines toward NGC 4945 is sig- 
nificantly larger than that in previous studies (Henkel et 
al. 1990; H94; COl). While the most complete study car- 
ried out so far, that of COl, contains a comparable num- 
ber of lines from CO and its isotopomers as well as from 
H2CO, the number of CS and SO lines studied in our sur- 
vey is larger. About 51% of the frequency range between 
81.4 GHz and 115.8 GHz, 42% of the frequency range be- 
tween 128.6 GHz and 159.5 GHz, and 27% of the fi-equency 
range between 217.7 GHz and 268.1GHz are covered. Of 
particular importance is HC3N which was not clearly de- 
tected in any previous study (see H94 for a tentative de- 
tection) but which was detected in our survey in at least 
half a dozen of lines. This molecule is useful to determine 
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Fig. 3. HNCO, HCN and HNC spectra measured toward the center of NGC4945 (see also the captions to Figs. 1 and 
2). 



H2 densities (e.g. Mauersberger et al. 1990) and can be 
used to check previously determined physical parameters 
of molecular clouds (e.g. COl). The (tentative) detection 
of a number of additional molecular species, like CH3CCH 
and CH3CN provides a broader view of chemical abun- 
dances while the detection of C^^S and the non-detection 
of -"^^CN give new insights into isotopic ratios in the nuclear 
environment of NGC4945 (see Henkel & Mauersberger 
1993 and COl for earlier summaries on extragalactic CNO 
isotope ratios, including NGC4945). 

4. Radiative transfer calculations 

With the peak and integrated line intensities determined 
from our measured profiles, we can derive a number of 
relevant line intensity ratios that can be used to estimate 
density, column density and kinetic temperature. The ex- 
citation problem involves statistical equilibrium of a mul- 
tilevel system that was solved in most cases by the Large 
Velocity Gradient (LVG) method. This allows us to deter- 
mine a local source function in which optical depths are 
related to escape probabilities of the photons. The model 
gas cloud has spherical geometry and uniform kinetic tem- 
perature and density. While there may be molecules which 



arc significantly affected by the intense radiation field in 
the nuclear region of NGC 4945 (for effects on HNCO, 
see Churchwell et al. 1986 and Sect. 4.9; for the negligible 
effects on CO, see Mao et al. 2000 and their Appendix 
C), such effects are difficult to quantify. If not said oth- 
erwise, a background temperature of 2.73 K is therefore 
used throughout the discussion. 

With lineshapes being approximately similar (for a 
discussion of the moderate differences obtained, see e.g. 
Bergman et al. 1992 and H94), measured peak intensi- 
ties (i.e. not the peak intensities obtained from Gaussian 
fits) appear to be more reliable input parameters for the 
model than integrated intensities (see Sect. 4.1). Taking 
CO source sizes according to Table 4 of COl, i.e. 29" for 
the J=l-0, 20" for the 2-1 and 15" for the 3-2 transi- 
tion, we corrected for beam dilution effects by calculating 
T4b=^mb/%f, with 7?bf = 0!/{0s+0l), and 9^ denoting 
the Gaussian FWHP source and beam size, respectively. 
Again following COl, for the other molecular species a 
source size of 20" was adopted. This appears to be a rea- 
sonable guess in view of the CO data and the extent of the 
1.3mm continuum emission (see Sect. 2.2). Table 4 gives 
beam dilution corrected line temperatures, which are used 
for the LVG calculations. 
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Fig. 4. HC3N, CN and OCS spectra measured toward the center of NGC 4945 (see also the captions to Figs. 1 and 2). 



Results of the model calculations are shown below, dis- 
playing lino intensities and line intensity ratios as a func- 
tion of H2 density and molecular column density. Common 
to all these plots are source size averaged line tempera- 
tures. Cloud averaged line temperatures may be signifi- 
cantly larger, depending on the degree of clumping. To 
give an example: if the HCN J=l-0 emission arises from 
the innermost 20" of the galaxy and if we assume a rota- 
tion temperature of rrot=10K and r (HCN l-0)> 1 (H94), 
the 'clumping factor' within the 20" sized region would 
become /c ~0.1. As long as the lines are optically thin, 
this does not imply a change in excitation or in the line 
temperature ratios. If, however, the transition from source 
to cloud averaged quantities involves optical depths near 
unity or higher, significant changes are expected. 

In the following, LVG simulations are discussed for 
some of the most important molecular species, carefully 
discriminating between optically thick and optically thin 
line emission. Most calculations were made for kinetic tem- 
peratures of 50 and 100 K. Such high values are justified as 
e.g. suggested by the multilevel CO studies of Mao et al. 
(2000) and Bradford et al. (2003) on M82 and NGC 253, 
respectively, and by the NH3 study of Mauersberger et 
al. (2003) towards NGC 253, Mafffei 2, IC342, and M82. 
As we shall see, even for Tkin = 100 K H2 densities are 



larger than previously anticipated (c.f. COl and Sect. 7) 
and would have to be further enhanced if lower Tkin val- 
ues would be adopted. 

4.1. CO 

Plotting CO/i^co J=l-0 and 2-1 line intensity ratios 
as a function of velocity (Fig. 6), we obtain maxima and 
minima that agree with those shown by Bergman et al. 
(1992; their Figs. 2 and 3) for the J=l-0 lines alone. This 
is consistent with the hypothesized ring-like morphology 
of the emitting gas (see also Whiteoak et al. 1990; Dahlem 
et al. 1993; Mauersberger et al. 1996a). 

Among the 10 detected lines of carbon monoxide (CO; 
Fig. 1, first three panels), the C^^O J=l-0 and the two 
J=3-2 profiles are the most uncertain. The signal-to-noise 
ratio of the C^^O J=l-0 line is small, but the line inten- 
sity is similar to that of the higher signal-to-noise detec- 
tion reported by COl. Our C^^'O J=2-l profile is of better 
quality. With respect to the CO and ^^CO J=3-2 profiles, 
that were not measured simultaneously, we note that the 
lineshape is narrow, mainly confined to the lower radial 
velocity range of the CO 1- and 2 1 profiles. Since we 
know from Mauersberger et al. (1996a) that CO J=3-2 
lines have no peculiar shape, the effect must be due to 
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Fig. 5. SO, CH3CCH, C2H, CH3CN, HCO+, 
also the captions to Figs. 1 and 2). 



SiO and HCO spectra measured toward the center of NGC 4945 (see 



small pointing errors. These can substantially affect our 
submillimeter data that were obtained with a relatively 
small beam. The CO and ^^CO J=3-2 profiles, while not 
showing emission over the entire velocity range, likely in- 
dicate the correct peak line temperatures, since the J=l-0 
and 2-1 line peaks are broad (> lOOkms"^) and since the 
peak temperatures arc not lower than those reported by 
Mauersberger et al. (1996a) and COl. This is the major 
motivation to use peak and not integrated line intensities 
for the model calculations. 

Applying an LVG code to the CO data (see Mao 
et al. 2000) and reproducing C^'^O and C^^O J=l- 
and 2-1 lines, we obtain an excitation parameter of 
Tkin • 't^hI^ ^ 2x 10^ K cm~^/^, i.e. a kinetic temperature of 
Ikin~100K yields a density of riHa ~3.5xl0'^cm~'^, in 
good agreement with COl (their Table 5). Also the result- 
ing column densities 

N = [X'/{dv/dr)] x 3.08 x 10^^ x Avi/2 

{N: column density in cm~^; X': LVG density pa- 
rameter of the specifically studied molecule in cm~^; 
dv/dr: velocity gradient in kms~^ pc~^; Avi/2- linewidth 



/2- 

in kms ^, the adopted value is 300 km s 



show 



with 3.7x1015 cm-2 and 2.4xl0i^cm-2 for C^^O and 
C^^O, respectively, good agreement. If we consider er- 
rors in individual line temperatures by ±10%, we obtain 
N{C^^O) /A^(Ci7O)=6.4±0.3. 

Fits to the ^^CO and CO peak line intensities, involv- 
ing three transitions in each case, are more difficult to in- 
terpret. Using the temperature and density deduced above 
implies 2 1/1 and 2-1/3-2 line intensity ratios that are 
larger than those observed. A ^'^CO column density of 
9.2 xlO^^ cm~^, that approximately reproduces observed 
line intensities, leads to 2-1/1-0 and 2-1/3-2 line inten- 
sity ratios of 1.8 and 1.6, respectively, that have to be 
compared with observed values of 1.3 and 0.9. Assuming 
that the 2-1 line, that has the highest optical depth, arises 
from a particularly large volume, makes the discrepancy 
even more pronounced (note that this is independent of 
the overall extent of the emission that is given in Table 4 
of COl; it refers instead to the source filling or clumping 
factor within that volume). The same holds for CO, where 
a model fit with N{CO) ^ LSxlO^^cm-^ yields 2-1/1-0 
and 2 1/3 2 ratios of 1.6 and 1.3, respectively, instead of 
the observed ratios of 1.2 and 1.0. Furthermore, an abun- 
dance ratio of ~20 between CO and ^^CO may be too 
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Table 3. A comparison of line intensities with COl, 
Mauersberger et al. (1996a) 



H94 and 



line 



J Tmb dv 

(K km s~^) 



Tmb 

(K) 



(2)" {Sr (4)" (5)" 



(7)'' 



CO 1-0 
CO 2-1 
CO 3-2 
i=*CO 1-0 
i=*CO 2-1 
i^'CO 3-2 
C^®0 1-0 
C^*0 2-1 
C"0 1-0 
C"0 2-1 
CS 2-1 
CS 3-2 
HCN 1-0 
HON 3-2 
HCO+ 1-0 
HCO+ 3-2 
H2CO 2i^2— 
H2CO 2o,2-lo,i 
H2CO 2i i-li 



523 
921 
871 
34.4 
81.2 
88.0 
9.29 
24.5 
1.21 
3.31 
17.0 
13.9 
22.4 
29.6 
24.0 
18.3 
3.73 
5.02 
5.78 



510 
740 
760 
30 
86 

8.4 
29 
1.8 

5 

9.3 
11 
24 
45 
21 
37 
6.5 
5.4 
5.9 



466'' 
1050 

32.3 
123.2 

10.9 
32.4 
1.33 



23.5 



1.8 

3.2 

3.5 

0.12 

0.24 

0.30 

0.03 

0.07 

0.004 

0.012 

0.04 

0.038 

0.065 

0.09 

0.08 

0.06 

0.013 

0.014 

0.019 



1.7" 

2.4 

1.6 

0.114 

0.37 

0.026 

0.12 

0.005 

0.019 

0.031 

0.038 

0.08 

0.2 

0.08 

0.225 

0.02 

0.017 

0.013 



ITS' 

3.2 

2.5' 

0.107 

0.4 

0.25'' 

0.037 

0.12 

0.0057 



0.08 



^ (1) lines selected from COl except "'^^CO 3—2, which is from 

Mauersberger et al. 1996a; 
'° (2) & (5) our data; 
■= (3) & (6) COl; 

(4) & (7) H94; 

^ applied beam efficiencies for columns (4), (6) and (7): 0.75 for 
HCN 1-0 and HCO+ 1-0; 0.70 for CO 1-0, ^^CO 1-0, C^*0 1-0, 
C^^'Ol-O and CS2-1; 0.60 for CS3-2 and the three H2CO lines; 
0.50 for CO 2-1, ^^00 2-1, C^*0 2-l, and C^''0 2-l; 0.40 for 
CO 3-2, ^^C0 3-2, HCN 3-2 and HCO+ 3-2; 

' Mauersberger et al. 1996a 



small in view of a ^^C/^^C ratio of ^-^50, that was sug- 
gested as a likely value by H94 and COl. 

What can resolve the obtained discrepancies? As al- 
ready mentioned, J~'i-2 peak line temperatures may be 
particularly uncertain. Here calibration errors can be more 
significant than for the other lines. Our CO J=3 2 peak 
intensity is larger by 30% than that reported by COl. 
Taking their spectrum would lead to a reasonable model 
fit. The shape of our ^"^CO J— 3-2 spectrum (Fig. 1) is 
unusual. It remains open whether the narrow feature 
at ~430 kms~^ (Tmb 0.4 K) or the wider component 
at higher velocities (Tmb 0.25 K) represents the proper 
peak intensity value. A good match with the model would 
imply Tmb ~ 0.12 K. Thus there could be a serious discrep- 
ancy between the model and the observational data that 
has to be addressed by future studies. 

Measured 2-1/1 line intensity ratios are more accu- 
rate. The CO J=2-l line of COl is weaker than the cor- 
responding line measured by us and would enhance the 
conflict between the model and observations. Having ob- 
served the line on many occasions, however, we believe 
that our calibration (and that of H94) is more appropri- 
ate. Otherwise, differences with previously reported J=l- 
and 2-1 intensities are negligible. 

One way to reconcile observed line intensities with 
those of the model is to assume clumping. As long as 
all line emission is optically thin, intensities will be pro- 
portional to /c X T, i.e. to the product of the clumping 
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Fig. 6. CO/^^CO intensity ratios as a function of radial 
velocity for the J=l-0 (upper panel) and J=2-l (lower 
panel) lines. 



factor (the line emitting area within the cloud size given 
in Sect. 4) and optical depth. For a given line intensity 
a decrease in the clumping factor /c would then imply a 
corresponding increase in optical depth. If r approaches 
unity, however, /c x t does not determine line intensi- 
ties any longer and each of the factors has to be ana- 
lyzed separately. The C^^O, C^^O, ^^^CO and CO column 
densities were derived assuming /c=l, reducing all opti- 
cal depths to their possible minimum. For CO, however, 
this approach is too simple. If fc is identical for all CO 
isotopomers and /c < 1, optical depths and column den- 
sities will be larger than the so far calculated values ig- 
noring clumping within the CO emission region. Higher 
opacities imply more photon trapping and thermalization, 
thus reducing differences between J=l-0, 2-1 and 3-2 
line intensities and leading to a better agreement with ob- 
served line temperature ratios. As already mentioned, ho- 
mogeneous cloud coverage yields A^(CO) ~ 2 x 10^* cm~^ 
(Fig.7a,b), iV(CO)/Af(i3CO)-20, and 2 1/1-0 and 2- 
1/3-2 line intensity ratios of 1.6 and 1.3. A cloud aver- 
aged column density of A^(CO) ~ 6x10^^ cm~^ implies 
7V(CO)/A^(i3CO) -50 (/c ~ 0.4; see Fig. 7c) and line in- 
tensity ratios of 1.4 and 1.2, that agree better with the 
observed ratios of 1.2 and 1.0, respectively. 
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While moderate (0.2< fc <1) clumping is thus likely, 
the procedure is not accurate enough to provide reliable 
values for f^, N{CO) and N{GO)/N{^^CO). In the fol- 
lowing subsections, we will therefore try to constrain the 
12q^13q isotope ratio in other ways. Assuming a constant 
carbon isotope ratio throughout the nuclear region, this 
wiU yield the proper iV(CO)/A^(^^CO) abundance and 
12(-;^13(~; isotope ratios. A comparison with the beam av- 
eraged CO column density will then also determine the 
clumping factor. 

4.2. CS 

Unlike CO with its high column density, optical depth 
and small critical density (where rates of collisional and 
radiative deexcitation become similar), carbon monosul- 
fide (CS) is characterized by smaller column densities and 
larger critical densities. While CO may arise predomi- 
nantly from warm layers at the surface of clouds that are 
irradiated by a strong UV radiation field (e.g. Mao et al. 
2000), CS emission likely originates from denser clumps 
of gas deeper inside the clouds where kinetic tempera- 
tures might be lower. While the study of COl was con- 
fined to the CS J=2-l and 3-2 lines, our LVG model fit 
(see Mauersberger & Henkel 1989) also includes the J=5- 
4 line that is crucial for an excitation analysis as well as 
the J=3-2 lines of two rare isotopomers, those of C^^S 
and 13 CS. 

At a kinetic temperature of Tkin'^lOOK, the 
J=2-l, 3-2, and 5-4 line temperatures lead to 
A^(CS)-4xlOi4cm-2 and nn^ 5x Iff* cm^^ (see Fig. 8). 
The CS column density is smaller, the H2 density is larger 
than reported by COl, by factors of five. Assuming a lower 
kinetic temperature does not significantly alter N{CS), 
but increases the density further. At rkin=25K, riHa would 
be in excess of 10^ cm^'^. The reason for the discrepancy 
between COl and our data lies in the fits to the J=2-l and 
J=5 4 hues. With the parameters suggested by COl, the 
J=2 1 line would be stronger and the J=5-4 line would 
be weaker than observed. 

The model calculations of COl and our results both 
indicate that CS lines, if arising from a homogeneous 20" 
sized region, are optically thin (r (CS 3-2) ~ 0.086). Our 
C34S and "CS spectra thus imply iV(CS)/iV(C34S) ~ 12 
and N{CS)/N{^^CS) ~ 50. If, however, clumping within 
the assumed 20" sized emission region is substantial and 
source filling factors become /c <0.1, CS line saturation 
may play a role, leading to higher N{CS)/N{C^^S) and 
A^(CS)/iV(i3cS) ratios. The latter also holds if our tenta- 
tive i^CS profile is an artefact. 
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Fig. 7. LVG calculations with the H2 density riRa (in 
cm~'^) as a function of column density N (in cm^^) for 
CO (dash-dotted lines) and its isotopomers C^^O (solid 
lines), C^^O (dashed Hues) and ^^CO (dotted lines) at 
rki„ = 100K and /c=1.0 (a), Tki„ = 50K and /c=1.0 (b) 
and Tkin = 100 K and /c=0.4 (c). /c is the clumping factor 
(see Sect. 4.1). Transitions and brightness temperatures 
are given for each of the 'vertical' linos. The 'horizontal' 
lines connect points with a specific line intensity ratio in 
the {N, riHa) plane. 



4.3. HC3N 

Cyanoacetylene (HC3N) is another probe of the higher 
density gas, providing a particularly large number of tran- 
sitions that are, unlike those of CO and possibly CS, al- 
most certainly all optically thin. We observed 9 rotational 



HC3N lines, 6 or 7 of which were detected (first two pan- 
els and the first spectrum of the third panel of Fig. 4). 
Including a total of 22 lines and assuming that the J=17- 
16 line was detected, the LVG fit (see Mauersberger et al. 
1990) gives for Tki„=100K ^■(HCsN) - 7xl0"cm-2 and 
riHa ^ 10^ cm~3 (Fig. 9), the latter in good agreement with 
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Fig. 8. LVG calculations for CS and its isotopoiners 
C^^S and 13CS at Tki„~100K (a), Tkin~50K (b) and 
Tkin- 25 K (c). 

the results from CS. For smaller kinetic temperatures, 
H2 densities would be even larger. The maximum opti- 
cal depth is T (HC3N 12-11) ~ 0.0024 (lOOK), or r (HC3N 
10-9) ~ 0.0030 (50K). The clumping factor would there- 
fore have to be small, of order fc ^ 0.005 or loss, to induce 
line saturation. According to the model, the J=21-20 line 
is an order of magnitude weaker than the tentatively de- 
tected J=17 16 feature. Therefore it is not surprising that 
the J=24-23 and 25-24 transitions remain undetected. If 
the tentative J=17-16 line profile is not an artefact, the 
145.6 GHz feature is the suporpostion of two lines with ap- 
proximately equal intensities, the J=16-15 line of HC3N 
and the 2o2~loi line of ortho-H2CO (Fig. 4). Otherwise, 
H2CO is dominating the line profile. 

4.4. SO 

Seven lines of sulfur monoxide (SO) were observed, two 
of which were detected (first two panels of Fig. 5). As a 
comparison with the HC3N J=10-9 and 12-11 line inten- 
sities shows (Fig. 4), the 100.1 GHz profile is dominated 
by HC3N 11-10, not by SO 45 - 44 emission. An LVG fit 
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Fig. 9. LVG calculations for HC3N at Tki„ ~ 100 K (a) and 
Tkin ~ 50 K (b). 



to the SO 32-2i and 43-32 lines (Fig. 10; for molecular 

constants and collision rates, see Tiemann 1974; Green 
1994) gives for Tkin=100K, 7V(S0) - 1 x lO^^^ cm-^ and 
riHa ^lO^cm"-^. For Tkin=50K, is about 50% larger 
while 7V(S0) is not significantly changed. For the two 
detected lines source averaged optical depths are of or- 
der 0.005. Densities are high enough to yield excitation 
temperatures that are approaching T^in- The computed 
brightness temperature of the ground state Iq-Oi transi- 
tion {v ~ 30 GHz) is ~ 35% of the lines detected by us. 
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Fig. 10. LVG calculations for SO at T^n ~ 50 K (sohd 
lines) and Tkin ~ 100 K (dashed lines). 
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4.5. HCN 

H94 presented an HCN J=l-0 map that showed that the 
emission is not spatially resolved by the ~50" beam of 
the SEST. COl detected the HCN J=l and 3-2 tran- 
sitions, while Chin et al. (1999) reported the detection of 
the J=l-0 lines of HCN, H^^CN, and HC^^^, indicating 
a surprisingly small ^'*N/"'^^N isotope ratio. Here we com- 
bine the three J=l-0 transitions with detections of the 
HCN and H^^CN J=3-2 lines and a highly tentative de- 
tection of the HCN J=:4-3 transition (Fig. 3). The unde- 
tected DCN J=2-l line is bracketed by two comparatively 
strong lines, those of CH3OH (3-2) and C^^S (3-2), so that 
even the baseline remains undetermined (see Fig. 3). 

LVG fits with collision rates from Green & Thaddeus 
(1974) are shown in Fig. 11. In M 82, the bulk of the HCN 
emission arises from gas that is not only denser but also 
cooler than the CO emitting gas (e.g. Mao et al. 2000). 
We therefore present data with Tkin=50 and 100 K to 
study differences in the model parameters. For rkin=10GK 
and assuming optically thin emission we obtain a den- 
sity of riHs ~ 1-5x10^ cm~^ and column densities of 
A^(HCN) - 2.6 X 10^'' cm-2, iV(Hi3CN) - 1.7x lO^^ cm'^ 
and iV(HCi5N)~9.3xl0i2cm-2. For rki„=50K we find 
instead ~ 3x 10^ cm~^. Excitation temperatures in the 
low- J lines are ^8 and 6K for Tkin=100 and 50 K, re- 
spectively. The lines are subthermally excited. Therefore 
the HCN J=4-3 line should be 3-4 times weaker than 
the J=3 2 line, which is consistent with the obser- 
vations. Optical depths are of order 0.05-0.1 for the 
main species. This implies that line saturation likely af- 
fects the HCN J=l line and possibly also the J=3 2 
transition. Assuming Local Thermodynamic Equilibrium 
(LTE) and taking excitation temperatures from the LVG 
simulations, iV(H^'^CN)^2.6xl0-^''cm~^ from the inten- 
sity of the J=l-0 line. With ^^C/^^Cr^bO (Sect. 6.1), 
7V(HCN)~1.3xl0i^cm-2. 

Apparently, the gas emitting at 710 km s^-'^ is more 
highly excited than that observed at lower velocities. In 
the J=3-2 line the HCN peak is observed at 710kms~^, 
while the feature is less dominant in the J=l-0 transition 
(see also Sects. 5.1.3 and 5.3). 

The HCN/H^^CN J=l-0 line intensity ratio of 20- 
25, reflecting an abundance ratio of ~ 15 in the optically 
thin limit, is similar to that found for the corresponding 
CO/^^CO line intensity ratios. This indicates similar op- 
tical depths in the CO and HCN main isotopomers, while 
the line emitting volume may be larger in the case of CO. 
The line intensity ratios between the main and the ^^C 
bearing species can be taken as lower limits to the ^^C/^^C 
isotope ratio. 




LOG (N(HCN)) 



Fig. 11. LVG calculations for HCN (solid lines) and its 

isotopomers H^^CN (dashed Hues) and HC^^N (dotted 
line) at Tki„~100K (a) and rkin~50K (b). 

3i.2^2i_i Hnes (for details of the code, see Henkel et al. 
1980), for rkin=100 and 50K. For Tki„=100K we ob- 
tain iV(ortho-H2CO) ~ lO" cm'^ and nn^ ^4x10^ cm'^. 
Optical depths are highest in the J=2-l lines but do not 
exceed 0.013. Excitation temperatures are below 40 K. 
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Fig. 12. LVG calculations for ortho-H2CO at Tkin ~ 100 K 
(solid lines) and Tkin~50K (dashed lines). 



4.6. H2CO 



Formaldehyde (H2CO) line profiles are displayed in the 
last two panels of Fig. 2. Fig. 12 shows results from our 
LVG calculations, including the observed Ka=l tran- 
sitions of ortho-H2CO, i.e. the 2i^2-li,i, 2i^i-li_o and 



4.7. OCS 

For carbonyl sulfide (OCS) we have a tentative de- 
tection of the J=9-8 line (Fig. 4). If the detection is 
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real, we obtain with the LVG code already used by 
Mauersberger et al. (1995) and assuming Tkin=100K 
7V(OCS)~5xl0i'*cni^2 (Fig. 13). As can be seen in 

Fig. 13, the column density is not significantly changed 
for Tkj„=50 and 25 K, while the H2 density is poorly con- 
strained. 



T^I^=100K 

T^,„ = 50K 

Tii„ = 25K 




12 14 16 

LOG (N{OCS)) 

Fig. 13. LVG calculations for OCS at Tkm^lOOK (solid 
line), rkin~50K (dashed line), and Tkin~25K (dotted 
line). 



7V(CH3OH)-5.5xl0i''cm-2 and nn^ ~ 1 x 10'* cni^^. For 
Tkin=25K, the column density is similar, but the density 
becomes a few times lO'^cm"'^. Optical depths are less 
than 0.03 for a source size of 20". 
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4.8. CAVgOAV 

Eight methanol (CH3OH) features were searched for, three 
of which were detected (first two panels of Fig. 2). The 
J==2-l profile at 96 GHz is composed of 4 transitions, and 
the J=3-2 feature at 145 GHz is composed of 8 lines. An 
LVG source code for methanol (CH3OH) was kindly pro- 
vided by S. Leurini. Einstein coefRcients were taken from 
Cragg et al. (1993), energy levels were adopted from Xu 
& Lovas (1997), and collision rates are from Pottage et 
al. (2001, 2002) for Tki„ - 20K and from D. Flower (priv. 
comm.) for higher Tkin. Because of the large linewidths ob- 
served towards NGC4945, CH3OH lines are often blended 
and cannot be observed individually. In our calculations 
we therefore added the line intensities of the four J—2- 
1 transitions near 96 GHz, of the eight J=3-2 transi- 
tions near 145 GHz and of the sixteen J=5-4 transitions 
near 241 GHz, respectively. The J=2-l and 3-2 lines are 
found within a velocity interval of ~80kms~^ (25 and 
40 MHz, respectively) that is small with respect to the 
linewidth of the source. The observed J=5-4 feature is 
only slightly wider than 300kms~* (Table 2). Thus only 
few of the 16 lines can contribute significantly to the 
total intensity. According to the model, the 5_i-4_iE 
transition should dominate. In view of the weak HNCO 
10o,io~9o,9 hnc (Fig. 3) we assume that the contribution 
from the HNCO llo,ii^lOo,io hue is negligible and that 
the observed profile is entirely due to CH3OH J=5-4 
emission. The results (Fig. 14) show that for Tkin=100K, 



Fig. 14. LVG calculations for CH3OH with Tki„-100K 
(a) and Tki„ ~25K (b). 



4.9. HNCO 

Isocyanic acid (HNCO), observed in six rotational tran- 
sitions (see Fig. 3), is a particularly interesting molecule 
because it can probe the far-infrared radiation field 
(A '^SSO/im) with its Ka=l ladder (see Churchwell et al. 
1986). As a consequence, the rotational population distri- 
bution of HNCO can be dominated by far-infrared radia- 
tion. 

For the radiative transfer calculations of HNCO (see 
Lapinov et al. 2001) collision rates were taken from S. 
Green ( |http: / /www.giss.nasa.gov/ data/mcrates/#hnco| | . 
Taking a background temperature of 2.73 K and ig- 
noring the llo,ii — 10o,io hne that is blended by a 
CH3OH feature (see Sect. 4.8), we obtain for Tkin^lOOK 
iV(HNCO)'-2.3xlO"cm-2 and nn^ 1.6x10^ cm-^ . 
Optical depths are less than 0.02. 

Fig. 3 shows the 4o,4-3o,3 (or 4o-3o) and 4i_4-3i,3 (or 
4_i-3_i) line profiles. While the former Ka=0 transi- 
tion is clearly detected, the latter Ka=l hne may be 
an artefact. The line temperature ratio is therefore ^5:1. 
Reproducing a line intensity ratio of 5:1 with the model, 
we find a background temperature of Tbg ~ 30 K (see the 
lower panel of Fig. 15), that can be taken as an upper limit 
to the actual background temperature in the inner 20" of 
NGC4945. Consistent with Churchwell et al. (1986), our 
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excitation temperatures then reach an equilibrium with 
the radiation field so that observed line intensity ratios 
fail to constrain H2 densities. The solution with Tbg=30K 
cannot provide more than a crude approximation. In con- 
trast to our solution with Ttg=2.73 K (Fig. 15a), lines with 
measured intensities do not intersect in Fig. 15b and are 
thus not properly reproducing our measurements. A re- 
alistic approach, however, also requires the introduction 
of a 'dilution factor' that remains observationally uncon- 
strained. 





8 
















/ / 

' 




^ 4-3=0.097 






(a) 




^5,-4-0.046 

6- 'S=Q.049 

7- 6=0.069^ 








/ ' /-V 




10-9=0.014^^— — 




6 






1 1-10=0.051 - 












4, 4-3, 3=0.020 X 








( [ / 

























V ..^"^ / 














6- 


5/5-4=1.065 




4 


6 


-5/4-3=0.505 ^v' -^^ijT-^^. . 








8 


1 1 


\ \ \ \ \ \ \ 


] 

\ 


H 1 1 1 1 1 1 — 










/ 4-3-0.097 






(b) 


1 ! 
^ \ % 


\ 


/ 5-4=0.046 

5-5=0.049 

7-6=0.069 

10-9=0.014 — — 




6 




\ \ 




11-10-0.054 












V '*i,t-3,.3=0.020 






























1 


1 


Tt,=30K 
Tkin=10CK 










13 14 15 16 17 

LOG {N(HNCO)) 



Fig. 15. LVG calculations for HNCO at Tkin ~ 100 K and 
Tbg~2.7K (a) and Tbg~30K (b). The latter case does 
not allow us to determine H2 densities, so that horizontal 
lines (near ~10^cm~^) are not shown. 



5. Other molecular species 

In Sect. 4 a variety of molecular species were analyzed with 
an LVG code. Other molecules, mostly with a smaller 
number of detected transitions, were analyzed with an 
LTE approach. In the following some of these molecules 
are introduced, among them CN with the N=l-Q and 2-1 
lines being observed. 

5.1. CN 

5.1.1. The [CN]/pCN] abundance ratio 

For the CN A''=l-0 transition wc find a group of ,7=1/2- 
1/2 lines (z/ ~ 113.17 GHz) and a group of J=3/2 l/21ines 
{v ~ 113.5 GHz) that were detected with high signal-to- 
noise ratios (Fig. 4). According to Henkel et al. (1998), CN 
provides the most stringent limits on ^^C/^^C isotope ra- 
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Table 4. Intensities used for LVG and LTE 
calculations 



Transition 


Tib 






(K) 


CO 1-0 


6.143 


0.293 


CO 2-1 


7.080 


0.452 


CO 3-2 


7.000 


0.500 


^^CO 1-0 


0.410 


0.293 


"CO 2-1 


0.531 


0.452 


"CO 3-2 


0.600 


0.500 


1 S /---V -I 

C O 1—0 


0.102 


0.293 


C^*0 2-1 


0.155 


0.452 


C^''0 1-0 


0.014 


0.293 


C^^O 2—1 


0.027 


0.452 


CS 2-1 


0.320 


0.125 


CS 3-2 


0.156 


0.243 


CS 5—4 


0.042 


0.471 


C S 3-2 


0.013 


0.237 


^^CS 3-2 


0.004 


0.222 


CH3OH 2-1 


0.180 


0.122 


CH3OH 3-2 


0.126 


0.238 


CH3OH 5-4 


0.054 


0.465 


H2CO 21.2-I1.1 


0.066 


0.228 


H2CO 2i i-li_o 


0.063 


0.252 


H2CO 3i.2— 2i,i 


0.014 


0.431 


HNCO 4o 4-3o,3 


0.097 


0.103 


HNCO 5o'5-4o'4 


0.046 


0.152 


HNCO eo^e-Sois 


0.049 


0.205 


HNCO 7o,7-6o,6 


0.069 


0.260 


HNCO 10o,io-9o.9 


0.014 


0.418 


HNCO llo^ii-lOo.io 


<0.054 


0.465 


HCN 1-0 


0.622 


0.105 


HCN 3-2 


0.176 


0.512 


HCN 4-3 


0.038 


0.651 


H"CN 1-0 


0.035 


0.100 


H"CN 3-2 


0.007 


0.499 


HC^^N 1-0 


0.020 


0.099 


HCsN 9-8 


0.088 


0.091 


HC3N 10-9 


0.055 


0.110 


HC3N 11-10 


0.062 


0.130 


HC3N 12-11 


0.080 


0.151 


TTPoN 1 ^ 14 
nv_'3 LO — J-^ 


037 


217 


HPqN 1 fi— 1 ^ 


058 


240 


HC3N 17-16 


0.030 


0.262 


SO 32-2i 


0.035 


0.128 


SO 4:5 32 


0.032 


0.221 


OCS !) >s 


U.(i:i:i 


(1.151 


CN 1-0 J=3/2-l/2 


0.436 


0.161 


CN 1-0 J=l/2-l/2 


0.263 


0.160 


CN 2-1 J=5/2-3/2 


0.162 


0.433 


HNC 1-0 


0.413 


0.109 


HN"C 1-0 


0.028 


0.101 


HCO+ 1-0 


0.797 


0.106 


HCO+ 3-2 


0.115 


0.520 



for the definition of ?7bfj see Sect. 4 



tios in cxtragalactic molecular clouds. Because the N=l-Q 
transition is split into several lines (spin-rotation and nu- 
clear magnetic and electric quadrupole interactions), opti- 
cal depths are not large and can be determined by an anal- 
ysis of line intensity ratios. Deviations from LTE (Local 
Thermodynamic Equilibrium) intensities appear not to be 
large. 

In order to determine the optical depth of the CN 
N=l-Q lines, we denote the J=3/2-l/2 lines with the 
index 1 and the J= 1/2-1/2 lines with the index 2. With 
an observed line intensity ratio of 1.58±0.06 (Table 2) and 
ti=2t2 (see Skatrud et al. 1983 for relative LTE intensities 
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under optically thin conditions) wc then find 



/ T^b,idv ^ 1 - e-^^ ^ 1 - e-"^ 



= 1.58 ±0.06. 



For n we obtain ti = 1.09±0.21 

The J=3/2-l/2 and 1/2-1/2 lines are part of the same 
spectrum so that their line intensity ratio is not affected 
by calibration errors. With ti and T2 known and 



and 



I2 oc 



TmhAdv- 



1 - e-"^! 



tmb,2 



dv ■ 



T2 



1 



where Sfifj,! and T,gi^2 are 0.667 and 0.333, respectively 
(Skatrud et al. 1983), wc then obtain for the 'opti- 
cal depth corrected' total CN iV=l-0 line intensities 
/i=J2=60.4±5.4 K km s-^ 

For -^^CN (spectra are shown in Fig. 4) we took the 
frequencies and relative line intensities (LTE, t<C1) from 
Bogey et al. (1984) and Gerin et al. (1984). With a 3a 
noise level of 5.7mK in an 18.82kms^^ wide channel (see 
Table 2 and Fig. 4) we obtain for a 300kms~^ wide line 
an upper limit to the integrated intensity of 



3a 



jTdv 



5.7 mK- 18.82 kms"^- ^300/18.82 



The resulting 3a upper limit is 0.43 K km s"-'^. 

As in the case of CN we have to correct for the fact 
that the ^^CN A''=l-0 lines are distributed over a fre- 
quency range (108.4-108.8 GHz). In contrast to CN, ^^CN 
does not require any corrections related to line saturation. 
The most prominent group of ^^CN lines, located between 
108.7802 and 108.7964 GHz (the difference corresponds to 
~45kms~^ that is small with respect to a total linewidth 
of 300kms~^) comprises 42.7% of the total LTE intensity. 
As a consequence we find 

[^2CN]/[^3CN] > 60.4/(0.43/0.427) ~ 60. 

Henkel et al. (H94) reported a 3a lower limit to the 
HCN/H^^cjs^ lij^g intensity ratio of 31, while the actual 
ratio is 23.4±1.9 (Table 2). In view of this discrepancy 
we conservatively estimate that the [-'^^CN]/[-'^^CN] abun- 
dance ratio should be ~45 or larger. 

5.1.2. CN line intensity anomalies 



The N=l-0 transition: 

As already mentioned, the CN transitions represent 
the superposition of a number of hyperfine (HF) com- 
ponents that are displayed in Fig. 16. While the N=l- 
J=3/2 1/2 profile ('HFl' in Fig. 16a) consists of a 
broad blue- and a narrow redshifted line, the J=l/2- 
1/2 profile ('HF2' in Fig. 16a) is wider and shows a single 
broad feature (see also the observed spectra in Fig. 4). Is 



this difference caused by a different intrinsic lineshape? 
A fit with two velocity components, characterized by 
i;i=450kms~\ 1=300 kms"^ and ?;2=710kms~\ 

^■(^1/2,2=50 km s~^ reproduces the observed J=3/2-l/2 
and 1 /2-1 /2 profiles in a satisfactory way (the linewidths 
refer to full width to half power Gaussians in optical 
depth). The lack of a prominent redshifted spike in the 
J= 1/2-1/2 lines is not caused by a different intrinsic line- 
shape but is a consequence of the different frequencies and 
strengths of the individual HF components (see Fig. 16). 

For the fit to the line profile we assumed optical 
depth ratios between individual HF components that are 
consistent with the relative intensities given by Skatrud et 
al. (1983). To scale optical depths, we took I]rj=1.09 for 
the J=3/2 1/2 and Sri=0.545 (ti=2t2, see Sect. 5.1.1) 
for the J=l/2 1/2 HF components, the index i denoting 
individual HF features. It is remarkable that the sum 
of these optical depths is, within the limits of noise, 
the same for the blue- (^450kms~^) and redshifted 
(~710kms~^) velocity components. 



The A''=2-l transition: 

A comparison of the iV=l-0 J=3/2-l/2 and iV=2- 
1 CN profiles was given by H94, who suggested that the 
relative strength of the A^=2-l 710kms~^ feature (see 
Fig. 4) is higher than that implied by the 7V=l-0 J—3/2- 
1/2 lineshape. Such an excess in the A''=2-l 710 km s"""^ 
feature is also indicated by our new data, making use of 
the VI and V2 components defined by the iV=l-0 line and 
subdividing the HF components into three groups, HFl, 
HF2, and HF3 (see Fig. 16b). The excess becomes apparent 
when comparing the observed line intensity ratio between 
the 710 and 450kms~^ components, ^1.9 (Fig. 4), with 
that of our computation, ~1.5 (Fig. 16). 

Concluding that the strongest CN A^=2-l component, 
the 710kms~^ feature, is enhanced, we find that line sat- 
uration cannot play a significant role in this transition. 
Otherwise, the feature would be weaker, not stronger than 
expected in our analysis assuming optically thin emission. 



5.1.3. CN column density 

A comparison of total line intensities gives ~40Kkms~^ 
for the CN iV=l-0 and ~30Kkms-i for the N=2-l line 
(see Table 2). Accounting for the optical depth in the 
iV=l-0 transition, the integral would be ~60Kkms~^ 
if we correct for saturation effects (Sect. 5.1.1), while 
line saturation may not be relevant for the 7V=2-1 line 
(Sect. 5.1.2). Correcting these integrated line intensities 
for beam size effects, i.e., dividing the integrated line 
intensities by i%i (see Sect. 4 and Table 4), we obtain 
~375Kkms~^ for the iV=l-0 line accounting for line sat- 
uration and '~70Kkms~^ for the N=2-l line. The line 
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Fig. 16. a) Upper panel: The computed CN A''=l-0 pro- 
file composed of the two velocity components VI and V2 
(see also Sect. 5.1.2). Observed CN profiles are displayed in 
Fig. 4. It is apparent that the strong spike in the J=3/2- 
1/2 line (HFl) at 710kms~^ has no similarly prominent 
counterpart in the lower frequency J= 1/2 1/2 line (HF2) , 
located at ^1500 1600 kms^^ in the rest frequency frame 
of the 3/2-1/2 (HFl) features, b) Lower panel: The CN 
N=2-l spectrum computed with the same velocity com- 
ponents VI and V2 that were used to simulate the 1- 
spectrum in the upper panel. For the 7V=2-1 transi- 
tion, three groups of hyperfine components (HFl, HF2 and 
HF3) had to be taken into account. Note that the com- 
puted intensity ratio of the peaks at 710 and 450 kms'^ is hiO(2-iv=o) I I}iCO{F 



from gas that is characterized by a recessional velocity of 
710kms~^. If this component has an A7^=2-l/l-0 ratio 
that is twice the average value (sec H94), this implies 
Tex ~4.1K (see also Sects. 4.5 and 5.3). 

5.2. HNC 

The J=l-0 HNC and HN^^C hues (Fig. 3) were detected 
with an intensity ratio of 11.4±1.4, i.e. well below the 
i^C/i^C lower limit derived by CN (Sect. 5.1.1). This im- 
plies that the HNC J=l-0 line must be optically thick. 
With a smaller ^"^C/^^C line intensity ratio than mea- 
sured for HCN (Sect. 4.5) and with only half the inte- 
grated intensity, HNC is likely arising from a more com- 
pact region than HCN and HCO"*". Assuming LTE and 
an excitation temperature of Trot ~8 K that was derived 
for HCN, a molecule with similar dipole moment and ro- 
tational constant, we obtain A^(HN^^C)~3.5x lO^^cm"^. 
With i2c/i3c~50 (Sect. 6.1), Ar(HNC)~1.8xl0i5cm-2. 

5.3. HCO^, HCO and thermal SiO 

While the J=l-0 profile of HCO"*" shows, as other lines, a 
dominant component at 450kms~^ and a satellite compo- 
nent at 710 km s"-'^, the relative strength of the features is 
inverted in the J=3-2 spectrum (see Fig. 5). Apparently, 
the gas seen at 710kms~^ is higher excited (see also 
Sects. 4.5 and 5.1.3). 

H13C0+ is blended with the HCO Iqi - Oqo 
and SiO J=2-l v=0 transitions (see Garcfa-Burillo 
et al. 2002). In order to provide meaningful con- 
straints on line intensity ratios, the three profiles 
shown in Fig. 5 (lower right corner) were simulated. 
Assuming that the 450 and 710 km s""*^ components 
of a given molecular transition have equal intensities, 
and linewidths of 200 and 150kms~^, respectively, we 
find that SiO might be strongest while HCO might 
be weakest. With /HCO(F=i-o)/-fHCO(F=2-i) ~0.3, 
for the HCO Iqi - Oqo J=3/2-l/2 lines we find 

^ 2.0 and 



qualitatively /hi3co+ /^hco(f=2-i) 



=2-1) 



'2.5. 



^1.5, while the observed spectrum (Fig. 4) shows a ratio 
^1.9. 



temperature ratio is 2-1/1-0^^0.185. With 



0.185 = 4 X e"^ X 



1 



1 



(e^" - i)-i - jc^y - i)-i 

(e^ - l)-i - (ev - l)-i 



a;=hz/io/kTex, i/io=113.386GHz and 2/=hzyio/k-2.73=1.99 
we then obtain for the excitation temperature 
Tex ~3.1K and for the column density, assuming 
LTE, Af(CN)~2.5xlO"cm-2. We note that the excita- 
tion temperature is very small, indicating densities well 
below 10^cm~3 (see e.g. Fuente et al. 1995). Most of the 
A^=l-0 710kms~^ feature is originating from the HFl 
group of HF components (see Fig. 2 and Sect. 5.1.2), i.e. 



The resulting HC0+/H13(;iq-f ij^e intensity ratio is 

41 and excitation temperatures become 4.7 K (450 km s"-'^) 
and 5.4 K (710kms~^ component) in the case of optically 
thin emission. We note that the decomposition of SiO, 
HCO and H-'^^CO"'" is too tentative to provide meaningful 
error estimates. A search for SiO J=3-2 v=0 emission 
that is less contaminated by other spectral features would 
be worthwhile. The HCO+/H-'^''CO+ line intensity ratio, 
being much more uncertain than other such ratios, is not 
included in Table 5. 

5.4. C3H2 

Profiles of cyclopropenylidene (C3H2) are shown in the 
right panel of Fig. 2. The 85 GHz 2i2 - loi Hne (Fig. 2) 
shows the double peak structure with one main component 
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at 450 and the other at TlOkms"^ that also character- 
izes a few other species, most notably CN (see Fig. 4 and 
Sect. 5.1). Higher excited C3H2 transitions (Fig. 2), the 
86 GHz 432 — 423 line and the superposition of the 150 GHz 
4o4 — 3i3 and 4i4 — 893 lines (separation:~60kms~^) show, 
however, an approximately Gaussian lineshape. With the 
partition function (e.g. http://spec.jpl.nasa.gov), the elec- 
tric dipole moment, asymmetric rotor line strengths and 
excitation above the ground state (e.g. Thaddeus et al. 
1985; Vrtilek et al. 1987; Lovas & Suenram 1989) known, 
a column density assuming LTE and optically thin emis- 
sion can be derived. 

Using the 2i,2 — li,o and 43,2 — 42,3 lines that both be- 
long to the ortho-C3H2 species and that were observed 
with similar beam sizes, an excitation temperature of 
23±5K is obtained. The total column density for a source 
size of 20" and an ortho-para ratio of 3:1 is estimated to 
be -3.3x101"* cm"2 

5.5. C2H 

While Table 2 and Fig. 5 show the parameters of a one 
component Gaussian fit to the 87 GHz iV=l-0 transition 
of C2H, the line is actually split into two components 
with an intensity ratio of —3:1. This is caused by hyper- 
fine sphtting (e.g. Tucker et al. 1974; Ziurys et al. 1982). 
Two groups of HF components are present, a redshifted 
one at ~87.32GHz and a blueshifted one at ~87.41GHz. 
Their LTE intensity ratio in the optically thin limit is 2:1. 
Our observed ratio is higher, likely not because of maser 
amplification of the stronger feature, but because some 
part of the weaker blueshifted HF component, that near 
710 km s"""^, is superposed onto the dominant more red- 
shifted HF component observed. A simulation of the entire 
profile as for CN (see Sect. 5.1.2) or for the blend of HCO+, 
HCO and SiO features (Sect. 5.3), assuming an intensity 
ratio of 2:1 for the group of HF components (Tucker et 
al. 1974; Ziurys et al. 1982), an intensity ratio of 1:1 be- 
tween the 450 and 710 km s"* features and linewidths of 
200 and 150 km s"*, respectively, provides a good fit. Our 
assumption that the emission is optically thin then yields, 
with Tex=20K, TVc^H ~ 6.8 x lO^^cm-^. 

5.6. CHsCCH and CH3CN 

Extragalactic CH3CCH and CH3CN was first reported by 
Mauersberger et al. (1991). While GH3CCH was detected 
in NGC 253 and M 82, CH3CN was only seen in NGC 253. 

In the case of NGC 4945 we also obtain a clear de- 
tection of CH3CCH, while CH3CN is not detected con- 
vincingly (Fig. 5). This may only be due to blending by 
the CS J=3~-2 line, i.e. a search for the J=5-4, 6-5 or 7~ 
6 transitions of CH3CN would be meaningful. Applying 
iV(CH3CCH)- 2.8 X 10i2cm-2(Kkms-i)-i • /cHsCCH 
(Mauersberger et al. 1991), we obtain a column density 
of 2.4x101'' cm-2. 



6. Isotope ratios 

6.1. ^^C/^^C and ^^0/^0 ratios 

For a complete evaluation of the column densities dis- 
cussed in Sects. 4 and 5, a knowledge of isotope ratios in 
the nuclear region of NGC 4945 is indispensable. Line in- 
tensity ratios of -'^^C and ^^G bearing isotopomers provide 
direct information on line saturation in the more abundant 
species if the ^^G/^^G ratio is known. Table 5 summarizes 
the observed ^^Cf^^C, ^^N/i^N, i^o/isq, i6o/17o and 
'^^S/'^'^S values based on integrated intensities. To inter- 
pret the ratios we assume in the following that they reflect 
the isotopic ratios in the case that the lines are optically 
thin or that they provide lower limits to the isotope ratios, 
if the line of the more abundant species is optically thick. 
Fractionation and isotope selective photodissociation may 
not be important in the nuclear region of NGC 4945 (H94) 
and are thus ignored. The data do not yet permit a sep- 
arate analysis of the 450 and 710 km s"* velocity compo- 
nents and are therefore based on peak temperatures and 
integrated intensities of the entire individual spectra. 

There is a total of 8 measures of the ^^C/^^C ratio 
(Table 5). Four of them yield small values (10-15), two 
intermediate values (15-30) and two high values (^45). 
The comparison of line intensity ratios suggests that GO 
and HNG rotational lines have highest optical depths, fol- 
lowed by those of HCN. While it may appear strange, 
at first sight, that r(GO J=2-l)~ r(HNC J=l-0), HNG 
is presumably arising from a smaller area, where average 
H2 column densities are larger than those in the region 
being responsible for the bulk of the GO emission. Note 
that GO/i'^CO line intensity ratios decrease with increas- 
ing rotational quantum number J. This is expected in a 
warm molecular medium where CO column densities in in- 
dividual states rise almost proportional to their statistical 
weights (2J-I-1). 

With ratios or lower limits of 40-50, the ^^C/^^G val- 
ues derived from CN and GS are closer to the actual 
12q^13q isotope ratio. While CN did allow us to estimate 
optical depths, CS may show a direct albeit low signal-to- 
noise detection of its ^^C bearing species. 

Following H94, we obtain for the CO J=l-0 lines 

i2c/13g< [/(CO)//(13cO)] r(CO) - (15±2)t(CO) 
and 

1^0/1^0 < [/(GO)//(Gi^O)] r(GO) - (58±9)r(GO), 

where the inequality on the left part of both equations 
accounts for the fact that CO, the main isotopomer, is 
optically thick. 

For the J=2-l lines we obtain correspondingly 

i2c/i3c< [/(C0)//(i3C0)] r(CO) - (12±2)t(CO) 

and 

i60/i»0 < [/(GO)//(Gi**0)] t(GO) - (41±7)r(GO). 
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Table 5. Intensity ratios of isotopomers 



Molecules 


Transition 


Ratiol*) 


Ratio2'') 


co/^-"co 


1-0 


15.2±0.5 


15.0±2.1 




2-1 


11.3±0.1 


13.3±1.9 




3-2 


9.9±1.1 


11.7±1.6 


CN/i3CN<« 




>45 






3-2 


49.0±17.1 


39.0±5.5 


HCN/H^^CN 


1-0 


23.4±1.9 


17.8±2.5 




3-2 


30.0±5.2 


25.1±3.6 


hnc/hn"c 


1-0 


11.4±1.4 


14.8±2.1 


HCN/HC^^N 


1-0 


48.3±5.1 


31.1±4.4 


CO/C^^O 


1-0 


56.3±2.5 


60.2±8.5 




2-1 


37.5±0.8 


45.7±6.5 


co/c"o 


1-0 


433.0±83.2 


438.8±62.1 




2-1 


278.3±26.5 


262.2±37.1 




1-0 


7.7±1.5 


7.3±1.0 




2-1 


7.4±0.7 


5.7±0.8 


CS/C^*S 


3-2 


10.4±0.8 


12.0±1.7 



^' Ratios of integrated intensities (see Table 2). 
Corrections for beam size effects are less than 10%. 
Ratios of corrected peak intensities divided by r/bf (see 

Sect. 4 and Table 4) 

If ^'^CS or H^'^CO^ is not detected, the ratios arc lower 
limits. For HCO"*" only tire 710 kms^^ component was 
taken to derive the ratios. 

See Sect. 5.1.1 



If both ^^CO and C^^O were optically thin, we could 
then derive, neglecting fractionation effects and isotope 
selective photoionization (sec H94), ^^c/i^c ^ (0.26 ± 
0.05) X (i^O/isO) for the J=l-0 lines and ^"^G/^^G ~ 
(0.29 ±0.07) X (160/180) for the J=2-l lines. If, however, 
also i"^CO is optically thick, the values would be lower 
and would be closer to i^c/i^C - 0.10 x (i^O/i^O) that 
is found in the nuclear region of the Milky Way. 

As already mentioned, low J CO column densities al- 
most rise with their statistical weight, 2J+1, in a warm 
approximately thermalized medium. As a consequence, 
t(13C0 J=2-1)> r("CO J=l-0). The fact that the cor- 
relation between i^C/^^C and ^^O/^^O ratios is almost 
independent of the choice of transition (J=l-0 or J=2- 
1) does not prove but strongly suggests that the ^'^CO 
opacities are small (see also Bergman et al. 1992 and 
H94). This implies that the ratio between i^C/^^C and 
i6Qy'i8Q ratios does indeed drastically differ from that in 
the central region of the Galaxy. Giving a higher weight to 
the result from the J=l-0 lines because of a lower degree 
of i"^CO line saturation and smaller errors, we then find 
1^0/130^40 and i^O/i^O^lSO. 

Inflowing gas through a bar (Ott et al. 2001) might 
provide fresh material with i^C/i'^O in excess of 50, if we 
use the analogy to the Galaxy with a i^C/i^C isotope ra- 
tio of ~70 near the solar circle (e.g. Wilson & Rood 1994). 
However, in the middle or at the end of a violent starburst 
(see the conflicting statements of Koornneef 1993 and Ott 
et al. 2001 on future star formation in NG0 4945), i^O 
production might be even more effective than i^O produc- 
tion (i.e. 1^0/1^0 ratios may decrease) if metallicities are 
solar or higher (Henkel & Mauersberger 1993; Henkel et al 
1994b; Langer & Henkel 1995). With a i^O/i^O ratio of 
^^250 in the central region of the Milky Way we then find 
i2C/i3C~40-60 and i6O/i®O~150-240 for NGC4945. A 



detection of i^CN in NG0 4945 could strengthen the case 
for our upper limits to the i^C/i^C and i^O/i^O ratios. 

6.2. ^^O/^^O, ^^N/^^N and^^S/^^S ratios 

An 1^0/1^0 ratio of 6.4±0.3 was already derived in 
Sect. 4.1. This is the best known isotope ratio in the nu- 
clear region of NGC4945. Its value is higher than that 
measured in the interstellar medium of the Milky Way or 
in the Sun (i*0/i ''0=3.5 and 5.5, respectively; see Penzias 
1981) and can be explained by the peculiar isotopic com- 
position of the ejecta from massive stars that play an im- 
portant role in a starburst environment. 

For the central region of NGC4945, Chin et al. (1999) 
determined the i^N/i'''N isotope ratio on the basis of three 
HON J=l-0 lines, those of HON, Hi^ON and HOi^N. 
Their measured Hi^CN/HCi^N line ratio of ~2 combined 
with a i2C/i-''C ratio of -50 resulted in i^N/i'^^N-lOO, a 
value that is much smaller than any value measured in the 
galactic interstellar medium. As a consequence, massive 
stars must be a source of i^N. With an Hi^ON/HOi^N line 
temperature ratio of 2.1±0.3 (see Table 2) and a i^O/i^O 
ratio of 50±10 we obtain i^N/i^N=105±26, in good agree- 
ment with Ohin et al. (1999). 

For the disk of the Galaxy (galactoccntric distances 
3kpc< Dec <9kpc) Chin et al. (1996) found ^^S/^^S ra- 
tios of (3.3±0.5)(DGc/kpc)-h4.1±3.1. The ^^S/^^S ratio 
in NGC4945 can be directly obtained from the CS/O^^S 
J=3-2 line intensity ratio, since CS/i^CS may be close 
to the i^O/i"^C isotope ratio so that OS line saturation 
should not be very significant. In the 'worst' case, i.e. if the 
CS/i^CS intensity ratio is 40 while i2C/i3C~60 (Table 5 
and Sect. 6.1), the CS/C"^^S line intensity ratio has to be 
multiplied by 1.5 to obtain ^^S/'^^'S. As a consequence we 
estimate 3^S/3^S=13.5±2.5. This corresponds to the value 
determined by Ohin et al. (1996) in the innermost galac- 
tic disk {Dqc ~3 kpc). As a result we may have found not 
only in the inner galactic disk (Chin et al. 1996) but also 
in a nuclear starburst environment a ^^S/^'^S ratio that 
is substantially smaller than that in the solar system. So 
far this is the only nuclear starburst that was studied in 
32g^34g ^Yic upper limit to the i^O/i^O ratio is not 
strongly constrained (Sect. 6.1). Therefore a clear state- 
ment on overproduction, compared to ^^S, by massive 
stars cannot yet be given and more observational data are 
needed to establish a clear picture. A summary of isotope 
ratios is given in Table 6. 

7. Column densities and densities 

In Sect. 4 column densities for several molecular species 
were determined. Making use of the isotope ratios de- 
termined in Sect. 6 and adding several molecular species, 
where measurements did not yet justify an LVG treat- 
ment (Sect. 5) provides a total of 24 molecular column 
densities adopting Tkin=100K and source sizes as out- 
lined in Sect. 4. The column densities do not change sig- 
nificantly for Tkin=50K. Changes in source size would. 
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Table 6. Isotope ratios' 



'> in the case that Une saturation effects play a role. While 

this is unlikely in the case of HC3N (Sect. 4.3) and SO 
'(Sect. 4.4), this is less clear for N2H+, HCO+, H2CO, 
HNCO, CH3OH, C3H2 and CH3CCH. Detections of rare 
isotopomers would be needed to clarify the situation. 

Comparing our colnmn densities with those of COl we 
note that some of the values in Table 7 are significantly 

Data taken from Penzias (1981), Henkel&Mauersberger (1993), HenkelSmaller than thoSC qUOtcd by COl. While in the Case of 

et al. (1993, 1998), Wilson & Rood (1994), Chin et al. (1996), Harrison^Q ^nd HCO+ differences are almost negligible, for CS, 

et al. (1999) and this paper. For the Large Magellanic Cloud (LMC), see ttI'^^a i re i i e 

Table 1 of Chin (1999). HON and H^'^CN the values din:er by about one order of 

Nuclear starburst environment 
Isotope ratios in the central kpc of the Galaxy 
'^^ Isotope ratios in the local interstellar medium 





NGC''' 
1!) lo 


NGC'"' 


M82''' 


Gal,<-' 
(.'vnlvs 


Solar'" 
Circle 


Solar 
Syst. 




50±10 


40±10 


>40 




77±7 


89 


16Qyl8Q 


195±45 


200±50 


>90 


~250 


560±25 


490 


18Q/17Q 


6.4±0.3 


6.5±1.0 


^8 


^-3.5 


3.6±0.2 


5.5 


"N/i^N 


105±25 




>100 


>600 


450±22 


270 


32g^34s 


13.5±2.5 








32±7 


22 



Table 7. Coin 



itios and donsiiicp'^^ 



Molecule 


N 






(cm- = ) 


(cm-3) 


CO 


4.6x10^" 


4.0x10^ 


i^'CO 


9.2x10^^ 


4.9x10=* 


ci»o 


2.4x10^" 


2.6x10=* 


C^^O 


3.7x10^^ 


4.0x10=* 


CN 


2.5x10" 


'^l.OxlO* 


CS "1' 


3.7x10" 


5.0x10* 


c^^s 


2.7x10^^ 


5.5x10* 


^^'cs 


^7.4x10^^ 


5.5x10* 


so 


l.OxlO" 


1.0x10^ 


C2Hf> 


6.8x10^^ 




HON 


1.3x10^^ 


1.5x10^ 


H^^CN 


2.6x10^=* 


1.5x10^ 


HC^^N 


9.3x10^^ 


1.5x10^ 


HCO+ 


^1.9x10^^ 




HNC 


1.8x10^^ 




HN^^C 


3.5x10^^ 




N2H+ 


4.2x10^^ 




OCS "> 


^5.2x10" 




ortho-HaCO 


1.0x10" 


4.0x10^ 


HNCO 


2.3x10" 


1.6x10* 


C3H2 


3.3x10" 




HC3N 


7.0x101=* 


1.0x10^ 


CH3OH ") 


5.5x10" 


1.1x10* 


CH3CCH '> 


2.4x10" 





For Tkin^lOOK (sec also Sect. 8); for as- 
sumed source sizes, sec Sect. 4. 

Assumed: ^-C/^'^C^SO (Sect. 6.1) 
'^^ Line saturation estimated by the relative 
strength of the fine structure components 
d) 32s/34s_;^3 5 (Sect. 6.2) 

Assumed: Optically thin emission (ques- 
tionable at least in the case of HCO+); 
for H2CO, COl estimated an ortho-para 
ratio ~2.3. 

Assumed: Tex=20K 

Assumed: Tex=5K (see Sect. 5.3) 
>>) Tex=8K, as for HCN 

See Sect. 5.6 



however, have a major impact with column densities in- 
creasing with decreasing extent of the molecular emission. 
The assumption of similar source sizes for all molecular 
species, close to those determined for CO and the 1.3 mm 
continuum (Sect. 2.2). is most reasonable; in view of miss- 
ing interferometric data. A summary of estimated column 
densities and densities is given in Table 7. 

While the column densities of ^^CS and OCS are based 
on tentative detections and are therefore upper limits, 
the column densities of several species are lower limits 



magnitude. We do not know the cause of this effect but 

note that excitation temperatures determined by us for 
HCN are small and that our CS and HCN data are more 
complete than those of COl. 

H2 densities are higher than expected. While densities 
of several 10"^ cm^^ deduced from CO and of '^10"' cm^"^ 
from CN and CH3OH agree with estimates given by H94 
and COl, most molecular species indicate riHa ~10^ cm~^. 
Densities for Tkin=50K instead of 100 K would be even 
larger. We thus conclude that the nuclear region of 
NGC 4945 contains a prominent high density component 
that is most convincingly supported by our multiline anal- 
ysis of HC3N (Sect. 4.3). While the uncertainty in T^in and 
the choice of spherical or plan-parallel cloud geometry may 
introduce deviations by factors of up to three from the 
values given in Table 7, we nevertheless find that most 
densities are significantly larger than those estimated by 
COl. 



8. Relative abundances 

8.1. iV(H2), M(H2) and X {CO) 

Since the column density of CO is relatively well deter- 
mined for the inner ~20" of NGC 4945, an H2 density can 
also be calculated adopting a [C0]/[H2] abundance ratio 
of ~ 8 X 10-5 (see Bradford et al. 2003). The total column 
density is then iV(H2) ~ GxlO^^cm-a 

To check the reliability of the H2 column density, we 
use our measurements of the 1.3 mm continuum. This 
continuum, described in Sect. 2.2, appears to arise from 
a similar volume than the molecular gas and is thus in- 
terpreted as emission from dust grains associated with 
molecular clouds in the nuclear environment. For an esti- 
mate of A''(H2) we have to account for the contamination 
of the measured 'continuum' by the CO J=2-l emission 
line (all other lines contribute much less). With a sensi- 
tivity of 25Jy/K, the 920Kkms-^ correspond to a flux 
of 17600 Jy MHz. Taking the width of the bolometer filter 
(30 GHz), the contamination of the measured continuum 
by CO ./ = 2 - 1 in the central beam is 0.60 Jy or 43% of 
the measured fiux. This is larger than the 14% determined 
for the center of M82 (Thuma et al. 2000). Assuming that 
a contamination of 43% is representative for the central 
emission in NGC 4945, the pure continuum flux in the cen- 
tral beam is 1.3 Jy and 1.6 Jy in the mapped region. 
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Following Mezger (1990), the H2 column density can 
be obtained from 



1 



iV(H2) = 0.97 X lO^^cm^^ , xo.^,., X, 

^ ' (6'b/arcsec)2(Z/Z0)6Td x 

Here, is the beam size. Si, denotes the flux density in- 
tegrated over the beam, Z represents the metallicity, 
is the dust temperature, and x = 14.4/AnimTd- For the 
nuclei of external galaxies, 6=1.9 (this dimensionless pa- 
rameter allows for changes in the dust absorption cross 
section in various environments; see e.g. Mezger et al. 
1987) may be justified, and should be of order 50 K 
(IRAS 1989). If we assume that the metallicity is twice 
the solar value (for galactic disk abundance gradients, see 
e.g. Hou et al. 2000), our pure continuum flux density 
in the central beam (1.4 Jyx0.57'~0.8 Jy) translates into 
A''(Il2) = 4.5 X 10^2 cm~2, in good agreement with the 
determination from CO. The main uncertainty, which we 
estimate to be a factor of two in both directions, are the 
metallicity, the factor b, and only to a lesser degree the 
dust temperature. 

If we assume a column density of A'^(H2) ~ 
6x10^2 cm~2 (this also compares well with X-ray data 
suggesting A'^(H) ~ lO^^cm"^ for the nuclear environ- 
ment, while iV(H) ~ 4xl02^cm~2 is found for the di- 
rect linc-of-sight to the nucleus (Done et al. 2003)), 
X{CO)=N{}i2)/jTcoi-odvri-j^ (sec Tabic 4) becomes 
-3.5xl0i^cm"2 [j^ km-^]-^ (cf. Maucrsberger ct al. 
1996a). This is a factor of ~6.5 smaller than the stan- 
dard galactic X(CO) factor of 2.3x1020 cm~2 [Kkm"!]"!. 
Such a low conversion factor is consistent with values 
found in the central region of other starburst galaxies (e.g. 
Mauersberger et al. 1996a, b; Solomon et al. 1997; Weifi 
et al. 2001b). 

Table 8 gives fractional abundances of molecular 
species with respect to A''(H2) studied in the starburst 
galaxies NGC 4945, NGC 253 and M 82 and toward proto- 
typical galactic molecular clouds, the Orion Hot Core, the 
Orion Ridge and TMC-1. 

8.2. NGC4945, Orion and TMC-1 

Agreement between the fractional abundances derived 
from the central region of NGC 4945 and the galactic 
sources is obtained for CS, N2H+, H2CO, HC3N, and 
CH3CCH. A multilevel HC3N study of the Sgr A molecu- 
lar clouds near the center of the Milky Way provides, on 
smaller linear scales than for NGC 4945, similar excitation 
conditions, i.e. riHs ^ 10^ cm~^ for Tkin=80K (Walmsley 
et al. 1986). For SO and HCN, the abundances agree with 
those for the Orion Ridge and TMC-1, but are lower than 
the enhanced abundances of the Orion Plateau and Hot 
Core, respectively. The C2H abundance seems to be much 
higher in NGC 4945 than in the galactic sources, but ex- 
citation and optical depths are poorly constrained. The 
HNCO abundance in NGC 4945 is consistent with those 
of galactic sources containing warm (Tkin > 10 K) gas, 
while the abundance of cyclic C3H2 seems to be similar to 



that obtained in the TMC-1. The CH3OH abundance is 
higher than in cold quiescent galactic clouds, but is lower 

than in the Orion compact ridge. 

With respect to warm galactic sources, HNC is sig- 
nificantly enhanced. In the galactic disk, interstellar 
A^(HNC)/iV(HCN) abundance ratios cover a range of at 
least two orders of magnitude. In quiescent dark clouds 
the abundance ratio is close to or larger than unity (e.g. 
Churchwell et al. 1984; Harju 1989), while in warmer and 
denser clouds associated with massive star formation the 
ratio can decrease to values of ~0.01 (e.g. Goldsmith et al. 
1981, 1986; Schilke et al. 1992). It is thus highly surprising 
to find iV(HNC)/iV(HCN)~l in the warm environment of 
a nuclear starburst, like that in NGC 4945. 

8.3. NGC 4945, NGC 253 and M82 

NGC 253, M82 (NGC 3034) and NGC 4945 are with 

D=2.5 4.0 Mpc the three nearest starburst galaxies. All of 
them show a rich molecular spectrum and infrared lumi- 
nosities of a few times 10^° Lq . Line intensities in NGC 253 
and NGC 4945 can be similar while those towards M 82 
tend to be slightly lower (compare e.g. our Table 2 with 
Table 1 in Sage & Ziurys 1995). While the starburst in 
NGC 253 is compact, mostly originating from a 60 pc sized 
region slightly southwest of the nucleus (Telesco et al. 
1993), the starburst in M82 is more extended exhibiting 
two lobes at a distance of ~150 250 pc from the dynam- 
ical center that may be part of a molecular ring (Mao et 
al. 2000). 

Not only with respect to size and overall morphol- 
ogy but also with respect to molecular abundances, the 
nuclear regions of NGC 253 and M82 are quite distinct. 
While abundances of CN, CS, HCN, HNC, C3H2 and 
CH3CCH are similar, SiO, HNCO, CH3OH and CH3CN 
are clearly underabundant in M82 (see Mauersberger et 
al. 1991, 2003; Nguyen-Q-Rieu et al. 1991, 1992; Henkel 
et al. 1993; Mauersberger & Henkel 1993; Hiittemeister 
et al. 1997; Garcia-Burillo et al. 2001). Recently, an un- 
derabundancc in NH3 was also demonstrated (Weifi et al. 
2001a; Mauersberger et al. 2003), while HCO was detected 
in M82, not in NGC 253 (Garcfa-Burillo et al. 2002). 
According to Takano et al. (1995), molecules formed by 
ion-molecule and neutral-neutral reactions with a low ac- 
tivation energy barrier are seen in both galaxies, while 
species requiring neutral- neutral reactions with significant 
activation energy or dust grain related reactions are abun- 
dant only in NGC 253. 

In order to explain the chemical discrepancies, Mao et 
al. (2000) proposed the presence of a warm diffuse inter- 
stellar medium that is responsible for the bulk of the CO 
emission in M 82. This medium is a consequence of cloud- 
cloud collisions, shocks, high gas pressure, high stellar den- 
sities and numerous young massive stars that may have led 
to the evaporation of vast amounts of dense cool gas cre- 
ating a highly fragmented molecular cloud debris. While 
the use of their LVG code does not necessarily lead to an 
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Table 8. Logarithmic fractional abundances'^''^ 



Molecule 


NGC 
4945 


NGC 
253 


M82 


Orion 
Hot 
Core 


Orion 
Ridge 


TMC-1 


CO"' 


-4.1 


-4.1 


-4.1 


-3.9 


-4.3 


-4.2 


^■^co 


-5.8 


-5.7 


< 5 7 










K A 


-D.4 












- / .Z 


- 1 .Z 










CN 


-8.4 


-8.6 


-8.4 


<-9 5 


-8.5 


-9.2 


cs 


-8.2 


-8.7 


-7.9 


-8.2 


-8.6 


-8.6 




-9.4 














<-Q Q 


-10.3 










SO 


-9.0 






^-7.7 


^-9.0 


-9.0 


C2H 


-6.9 




-6.8 


g-9.5 


-8.3 


-8.6 


HON 


-7.7 


-8.0 


-8.0 


-6.5 


-8.3 


-8.4 


H^^CN 


-9.4 


-9.6 










HC^^N 


-9.8 












HCO+ 


^-9.5 


-7.7 


-8.3 


-8.6 


-8.6 


-8.1 


HNC 


-7.5 


-8.3 


-8.7 


-9.0 


-9.3 


-7.9 


hn"c 


-9.3 


-().() 










N2H+ 


-10.2 


-9.9 


-9.8 




-10.2 


-10.0 


OCS 


^-8.1 


-7.5 


<-7.9 




-8.5") 




HaCO"' 


-8.7 


-7.9 


-7.5 


-7.6 


-7.8''' 


-7.9 


HNCO 


-8.6 


-8.4 


-8.1 


-8.2 


^-8.7 


-10.1 


C3H2 


-8.3 


-9.2 


-8.7 




-10.1 


-8.0 


HC3N 


-9.0 


-7.9 


-8.0 


-8.8 


-9.9 


-8.9 


CH3OH 


-8.1 


-7.9 


<-8.7 


-6.5 


-6.9'^' 


-9.1 


CH3CCH 


-8.4 


-7.8 


-7.3 




-8.5 


-8.7 



^' Most fractional abundances of the TMC-1 were taken from 
Pratap et al. (1997, their weighted mean; see also Leung ct al. 
1984; Millar & Freeman 1984; Madden et al. 1989). Orion data 
were taken from Blake et al. (1987) and Comito (2003). Fractional 
abundances of NGC 253 and M 82 (Af(H2) = 4.0x10^^ cm"^ 
and l.Ox lO^^cm"^, respectively; Mauersberger et al. 2003) 
are, whenever possible, based on spectra from the 30-m tele- 
scope on Pico Veleta/Spain (see Henkel et al. 1988, 1993, 1998; 
Mauersberger &; Henkel 1989; Mauersberger et al. 1990, 1991, 
1995; Nguycn-Q-Rieu ct al. 1991, 1992; Hiittomeister et al. 
1995, 1997; Harrison et al. 1999; Oikc ot al. 2004), comple- 
mented by higher excited lines measured with similar beam size 
at the HHT (Heinrich Hertz Telescope) on Mt. Graham/USA 
(Mauersberger et al. 2003). For NGC 4945, sec this paper (Af(H2) 
= 6x10^^ cm"^); for N2H+, see Sage & Ziurys (1995). 

^' In most cases excitation conditions were assumed to be the same 
in NGC 4945. XGC 253 and M82. This is a simplification in view 
of the discussion about HC3N in Sect. 8.3. The CN opacity cor- 
rection derived in Sect. 5.1.1 for NGC4945, -t-0.2dex, was how- 
ever not applied to NGC 253 and M 82 (for the CN data, see 
Henkel et al. 1988). While HCO'*' emission was assumed to be 
optically thin in NGC 4945 (Sect. 5.3), a small HCO+/h1='CO+ 
line ratio (-^10) and ^^C/^^C ~ 40 (Henkel et al. 1993) yield 
'^HOO+ 4 NGC 253. Thus for NGC 253 0.6 dex were added 
to the column density calculated assuming optically thin emission. 
Fractional NaH^ abundances were obtained from 12-m Kitt Peak 
data (Sage &c Ziurys 1995) with a large beam size, that reduces 
their beam averaged column density significantly. This is com- 
pensated by the high rotational temperature (Trot — 30 K versus 
8K in Table 7). Af(N2H+) in NGC 253 and M82 and iV(HCN), 
JV(HNC) and Af(C2H) in M82 may suffer from optical depth and 
excitation effects that cannot be quantified. For details concern- 
ing column densities in NGC 4945, see Table 7. From the OCS 
column densities, derived by Mauersberger et al. (1985) for 11" 
sized regions, 0.6 dex were suljiracted to adjust A^(OCS) to a 
source size of 20". While there exists an HC3N multilevel study 
for NGC 253 (Mauersberger ct al. 1990) similar to the one out- 
lined in Sect. 4.3 for NGC 4945, only the HC3N J-10-9 line was 
reported from M 82 (Henkel ct al. 1988). Adopting the excitation 
conditions of NGC 4945 yields a logarithmic fractional abundance 
of -8.7, while excitation conditions similar to NGC 253 yield —7.3 
for M 82. The mean of these two values is given. 
For the extragalactic sources, the fractional CO abundance was 
taken from Frerking et al. (1982). The corresponding abundances 
for the galactic sources are based on direct measurements (e.g. 
Blake et al. 1987). 
Orion Compact Ridge 

^' Galactic abundances include para- and ortho-H2CO. For the 
extragalactic sources ortho-H2CO column densities (Table 7; 
Hiittemeister et al. 1997) were upgraded by 0.1 dex. 



area filling factor that is similar to the volume filling factor 
and while not all of their PDR model input parameters are 
self-consistent, more recently published data convincingly 
confirm their proposed scenario. CO is highly excited and 
must arise predominantly from gas with kinetic tempera- 
tures well in excess of 50 K (Mao ct al. 2000). The rotation 
temperature of NH3 is, however, only ~30K (Weifi et al. 
2001a; Mauersberger et al. 2003), implying Tki„~50K, and 
HCN data suggest temperatures of order Tkin ^20-60 K 
(Seaquist et al. 1998). Weifi et al. (2001a) also find that 
the fractional NH3 abundance decreases from low values in 
the lobes to even smaller ones near the dynamical center. 

Overall 60 /xm/lOO i^m dust temperatures for NGC 253 
and M 82 are (as the total infrared luminosities) quite sim- 
ilar. They are of order 50 K for a modified Planck function 
Si, = e,yB,y{Tfi), where is the dust temperature, i?i/(Td) 
the Planck function and ei, the emissivity of the dust grain 
population at frequncy i^. e,^ oc ly^ , /3=1. In M82 it thus 
appears that NH3, HCN and other high density tracers 
arise from dense gas with kinetic temperatures close to 
the dust temperature, that is too low for the synthesis of 
significant amounts of SiO and CH3CN or for NH3 and 
CH3OH evaporation from dust grains. The lower density 
tracer CO, however, predominantly originates from the 
warm UV irradiated surfaces of molecular clouds. This 
view is further supported by the comparatively larger 
amount of molecular gas at Tkin ^150 K, about ten times 
as much as in NGC 253 (Rigopoulou et al. 2002). At these 
temperatures, NH3 is evaporated from dust grains. Since 
we see NH3 at such temperatures in NGC 253 but not 
in M82, the warm (Tkin ~150K) molecular gas in M82 
must be tenuous, i.e. column densities of individual molec- 
ular filaments arc not high enough to shield the molecules 
from rapid photodissociation. Since HCO is observed in 
the Milky Way at the interfaces between molecular and 
ionized gas, the non-detection of HCO in NGC 253 and 
detailed maps of HCO in M 82 (Garcfa-Burillo ct al. 2002) 
further strengthen the idea that the molecular region of 
M 82 can be viewed as a giant PDR. Apparently, the star- 
burst in M 82 represents a later stage of evolution than 
the starburst in NGC 253. 

The crucial question then is whether the nuclear region 
of NGC 4945 resembles more closely that of NGC 253, that 
of M 82 or whether NGC 4945 represents the first known 
example of a starburst in yet another evolutionary stage. 
Based on the weakness of B7 emission Koornneef (1993) 
coined the term 'Postburst Infrared Galaxy' for objects 
like NGC 4945. Therefore similar physical and chemical 
conditions might be expected in the nuclear environment 
of M82 and NGC 4945. Ott et al. (2001), however, identi- 
fied a bar-like feature in NGC 4945 and claimed that large 
amounts of molecular material can still reach the nuclear 
region providing sufficient fuel for an ongoing starburst. 

Table 8 provides some clues to the stage of the star- 
burst in NGC 4945. Apparently, fractional abundances of 
CH3OH and HNCO are close to those of NGC 253, i.e. 
NGC 4945 does not show the underabundances that char- 
acterize the molecular spectrum of M82. For NGC 4945 
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and NGC253, carbon and oxygen isotope ratios appear 
to be identical within the error limits, while the ratios 
are less well constrained for M82 (Table 6). An important 
test, the measurement of thermal SiO, HCO and CH3CN, 
is, however, still required to obtain a complete picture. 
While wc may have seen some traces of SiO in NGC 4945 
(see Fig. 5 and Sect. 5.3), better data are needed for a com- 
parison. 

An evaluation of molecular line excitation can be made 
for HC3N, comparing our data with the HC3N multilevel 
study of Mauersberger et al. (1990) towards NGC 253. In 
the latter galaxy, excitation is significantly higher than 
in NGC 4945, i.e. there is a high excitation component 
leading to a detection of the J=26-25 line. For NGC 4945, 
our model calculations (Sect. 4.3) clearly indicate that the 
line is undetectable even if a 30-m sized telescope had been 
used. 

8.4. CN, HNC and HCN in NGC 4945 and other 
galaxies 

Huttemeister et al. (1995) and Aalto et al. (2002) com- 
pared CN iV=l-0 and 2-1, HNC J=l-0 and HCN J=l-0 
line intensities in a representative sample of galaxies, in- 
cluding targets that resemble the Galaxy (e.g. IC342), 
weak starburst galaxies (e.g. NGC 253 and M82), lumi- 
nous infrared galaxies (LIRGs; e.g. Arp299) and ultralu- 
minous infrared galaxies (ULIRGs; e.g. Arp220). As al- 
ready mentioned in Sect. 8.2, galactic HNC/HCN abun- 
dance ratios vary by at least two orders of magnitude 
(Schilkc et al. 1992). This is, however, not fully seen in the 
extragalactic line intensity ratios, (1) because particularly 
high HCN abundances lead to optically thick J=l-0 lines 
with intensities that do not fully reflect the column den- 
sity and (2) because the surveyed regions of ^1 kpc in size 
smear out any extreme values on small spatial scales. Thus 
[HNC]/[HCN] < /(HNC)//(HCN) <, 1 with [X] and /(X) 
denoting the abundance (relative to H2) and intensity of 
the respective molecule. Equal excitation is assumed and 
is supported by similar molecular constants and electric 
dipole moments. 

Towards NGC 4945, the HNC abundance is high. 
/(HNC)//(HCN)-0.5 is near the average of the values 
determined by Hiittemeister et al. (1995) and Aalto et al. 
(2002). The ratio of the optically thin lines of the rare 
isotopomcrs HN^'^C and H^'^CN indicates, however, that 
[HNC]~10~^'^. This is exemplified by the relatively small 
HNC/H^^NC line intensity ratio in NGC 4945 (see Table 
2) and the relatively large one in NGC 253 (Henkel et al. 
1993), indicating lower optical depths in the latter case, 
since the ^^C/^^C carbon isotope ratios appear to be sim- 
ilar in both galaxies. 

As already mentioned in Sect. 8.2, high HNC/HCN 
line intensity and abundance ratios (near unity) 
are a characteristic property of dark clouds with 
small temperature (Tkin~10K) and moderate den- 
sity (n(H2)~10*~^cm~^). At higher Tkin values and 



Table 9. Relevant line intensity ratios 

CN Ar=l-0/HCN ,7=1-0 1.79 

CN Ar=2-1/HCN ,7=1 1.31 

HCO+ .7=l-0/HCN .7=1-0 1.07 

HNC J=l-0/HCN J=l-0 0.54 

HN"C jV=l-0/H"CN J=l-0 1.10 

Taken from Table 2. Uncertainties amount to ~±15%. 



densities HNC may more efficiently react e.g. via 
HNC-hO-*CO-|-NH and HNC+H-»HCN4-H that shifts 
the chemical balance between HNC and HCN in favor of 
the latter species. Rcicciiitly, however, Herpin & Cernicharo 
(2000) showed that in the warm (Tkin~250K) and dense 
(n(H2)^10'') circumstellar torus of the protoplanetary 
nebula CRL618 HNC and HCN are similarly abundant. 
Aalto et al. (2002) proposed that [HNC]/[HCN]~1 is 
not only typical for dark clouds but may also character- 
ize^ warm regions as long as the ion abundance is high. 
Protonized HCN and HNC, i.e. HCNH+, will then recom- 
bine to produce HCN and HNC in equal quantities. 

A cool molecular halo of moderate density with par- 
ticularly high HNC abundance that surrounds the nuclear 
region was alternatively proposed by Huttemeister et al. 
(1995). This is, however, not a viable explanation for all 
sources, since our study of NGC 4945 is confined to the 
inner ~20" and a cool moderately dense foreground layer 
from the outer disk (NGC 4945 is viewed almost edge-on, 
see Sect. 1) would have been seen in other molecular lines, 
most notably those of CO. Another alternative, a cool 
moderately dense gas component within the nuclear envi- 
ronment does also not explain all observations. While such 
gas is common in the central kiloparsec of our Galaxy (e.g. 
Hiittemeister et al. 1997), ULIRGs like Mrk231 do not ex- 
hibit a cool dust component that would be associated with 
such a hypothesized cool source of HNC emission. 

Aalto et al. (2002) outlined an evolutionary scenario 
starting with bright HCN emission from dense quiescent 
warm gas with little CN and HNC, reaching a PDR- 
phase with /(CN)>/(HCN) and ending with bright HNC 
but little CN, when the starburst has evolved beyond a 
strong radiative impact from young massive stars, while 
the cloud's chemistry is still dominated by ion-neutral re- 
actions at moderate gas density and relatively high elec- 
tron abundance. NGC 4945 is in the intermediate phase, 
with /(CN)>/(HCN) (see Table 9), like NGC 253 (see 
Henkel et al. 1993), but with considerably more HNC, 
suggesting a later stage of evolution. M82 is also char- 
acterized by /(CN)>/(HCN), while the optical depth of 
HNC remains poorly constrained (Henkel et al. 1998). The 
starburst in M 82 is in a later stage of evolution, because 
here it is not apparent that large amounts of fresh molec- 
ular fuel continue to reach the central region. The /(CN 
N=l-0 J=3/2-l/2)//(HCN J=l-0) ratio may serve as 
an indicator for the importance of PDRs (e.g. Greaves & 
Church 1996; Rodriguez-Franco et al. 1998). In the galac- 
tic nuclear region, the ratio is ^0.3 (Paglione et al. 2002). 



M. Wang et al.: Dense gas in nearby galaxies 



23 



For NGC253, NGC4945 and M82, we find 0.8, 1.1 and 
1.3, respectively. 

9. Energetics 

The heating of gas in the nuclear environment of spiral 
galaxies is one of the major puzzles in our understanding 
of the interstellar medium. A variety of potential heating 
mechanisms has been proposed but to discriminate be- 
tween the various physical processes has turned out to be 
extremely difficult. Proposed were (1) cosmic ray heating 
(e.g. Glisten et al. 1981; Bradford et al. 2003), (2) stellar 
ultraviolet heating via PDRs (e.g. HoUenbach & Tielens 
1997; Mao et al. 2000; Schulz et al. 2001; Garcfa-Burillo et 
al. 2002), (3) heating by X-rays, e.g. from the active galac- 
tic nucleus (AGN) resulting in an X-ray dominated region 
(XDR; e.g. Lepp & Dalgarno 1996; Usero et al. 2004), (4) 
dynamical heating by shocks (e.g. Martin-Pintado et al. 
1997, 2001; Hiittemeister et al. 1998; Garcfa-Burillo et al. 
2000, 2001) triggered by outflowing material from massive 
stars or active galactic nuclei, by cloud-cloud collisions, 
by tidal forces or by clouds gravitationally dominated by 
the mass distribution of the embedded stars (e.g. Downes 
et al. 1993; Mauersberger et al. 1996b) as well as (5) by 
ambipolar diffusion or ion slip that might be particularly 
signiflcant in the potentially strong magnetic fields of the 
nuclear environment of spiral galaxies (e.g. Scalo 1977). 

In local dark clouds, cosmic rays are a major source 
of energy, yielding kinetic temperatures of 5-10 K alone. 
With a detailed chemical model, also involving cooling by 
the fine structure lines of O and C and by CO, Farquhar 
et al. (1994) find for cosmic ray heating Tkin=6.7 S°-^^ K, 
with the factor S denoting the enhancement {S >1) of 
the cosmic ray flux over the solar system value. Adopting 
S'^IO^ for sources hkc NGC 253, M 82 and NGC 4945 (e.g. 
Bradford et al. 2003; for supernova rates, sec Ulvcstad & 
Antonucci 1997), we thus obtain Tkin^SO K, in good agree- 
ment with the temperature of the dust in the starburst 
environment of all three galaxies. This is consistent with 
a scenario for M 82 that assumes a significant contribution 
to the overall heating by cosmic rays for the dense deeply 
embedded gas, while the warmer CO is heated by PDRs 
(see Sect. 8.2). In NGC 253, however, both CO (Bradford 
et al. 2003) and NH^ (Mauersberger et al. 2003) arise from 
regions characterized by Tkin ^100 K, so that at least for 
these two molecules, cosmic ray heating appears not to 
be sufficient. So far, NGC 4945 is less thoroughly studied. 
Detections of NH3 and of highly excited CO lines have not 
yet been obtained. To confirm or reject our finding that 
NGC 4945 appears to resemble more NGC 253 than M 82, 
higher excited (J >3) CO lines and NH3 should also be 
measured. 

Since high J CO and NH3 line data are not yet avail- 
able, indications for high molecular cloud temperatures 
are still rare for NGC 4945. Perhaps the strongest hint for 
the presence of warm molecular gas at mm-wavelengths 
is our HC3N multilevel study (Sect. 4.3), where the as- 
sumption of a 'normal' density also seen in most other 



molecular tracers, riHa ^10^ cm""^, leads to a kinetic tem- 
perature of Tkin ^100 K. While the presence of gas at such 
temperatures is a typical property of PDRs, HC3N is a 
tracer of deeply embedded molecular material that may 
be less affected by UV heating than CO. Furthermore, 
HCO, a tracer of PDRs, appears to be weak in NGC 4945 
(Sect. 5.3). While PDRs will certainly exist in the nu- 
clear environment of NGC 4945, they appear not to be 
widespread enough to cause heating on large spatial scales. 
Heating by shocks or ambipolar diffusion are thus more 
likely the physical processes that may heat the gas on large 
scales in the nuclear starburst environment of NGC 4945. 
In view of the presence of shocks that are expected to be 
associated with the bar identified by Ott et al. (2001), the 
measurement of a well calibrated SiO J=3-2 {v=0) profile 
would be highly desirable. 

An alternative heating source could be X-ray radia- 
tion, leading to a high ionization rate (Lepp & Dalgarno 
1996; Usero et al. 2004). Here we expect that the 20" re- 
gion, corresponding to a linear scale of 390 pc, is too large 
to be dominated by the AGN, which is consistent with 
an /(HCN)//(CO) line intensity ratio that is far from 
the extreme value measured toward the central region of 
NGC 1068 (Sternberg et al. 1994). The observed abun- 
dances of CO, HCO+ and HCN match those calculated 
by Lepp & Dalgarno (1996) for an ionization rate per hy- 
drogen atom, ^/n(H), of ~10~^^ "^^® cm^'' s^^, but the 
measured CN abundance (see Table 8) is an order of mag- 
nitude too low. 

10. Conclusions and prospects for future research 

Toward the nucleus of the southern starburst galaxy 
NGC 4945, 80 lines from a total of 19 molecules have been 
observed, including 49 detected line features, 9 tentatively 
detected transitions and 22 undetected lines. 

Two velocity components, one at ~450kms~^, the 
other at ~710kms~^, characterize many but not all spec- 
tra. The former component is wider; the latter component 
is higher excited as is indicated by HCN, HCO+ and CN 
spectra. While no maps are presented, CO/"^^CO line in- 
tensity ratios as a function of velocity are consistent with 
a ring- like morphology of the emitting gas. 

Using the LVG or LTE approximation, H2 densities 
and fractional abundances of 24 molecular species (includ- 
ing isotopomers) were calculated. Many of these species 
indicate the presence of a prominent high density inter- 
stellar gas component characterized by nH2~10^cm~^. 
Excitation temperatures are in some cases quite low and 
reach 3-4 K in the case of CN. HNCO data indicate a 
mm- wave background radiation field of ^30 K. 

Based on column densities derived from CO and from 
the 1.3 mm dust continuum (2.8±0.3Jy), the H2 to in- 
tegrated CO J=l-0 line intensity conversion factor 'X' 
becomes ~3.5xl0^^cm~^ [K km~^]~^. This is a factor 
~7 smaller than values found for the galactic disk. 

Calculated fractional abundances are compared with 
abundances observed toward the starburst galaxies 



24 



M. Wang et al.: Dense gas in nearby galaxies 



NGC253 and M82 and selected galactic sources. The 
chemical properties of the gas in NGC 4945 resemble more 
those in NGC 253 than in M 82. There are, however, differ- 
ences between NGC 4945 and NGC 253. One is excitation. 
HC3N is more highly excited in NGC 253. Furthermore 
there is an 'overabundance' of HNC in the nuclear environ- 
ment of NGC 4945. The HNC/HN"C J =1-0 line intensity 
ratio is with '^11 particularly small indicating a high op- 
tical depth (~5) in HNC J=l-0. Thus NGC 4945 is one 
of the few known starburst galaxies with an HNC/HCN 
abundance ratio ~1 that may be caused by a high abun- 
dance of ions. A detailed analysis of CN, HCN and HNC 
line intensities and abundances suggests that the starburst 
in NGC 4945 is in a stage of evolution that is intermediate 
between those of NGC 253 and M 82. 

Carbon, nitrogen, oxygen and sulfur isotope ratios are 
determined. Carbon and oxygen isotope ratios in the nu- 
clear environment of NGC 4945 and NGC 253 appear to be 
almost identical, while ratios from M 82 are less well con- 
strained. High 1^0/1^0, low i^o/iSQ and i^N/i^N and 
perhaps also low ^^S/^^S ratios (6.4±0.3, 195±45, 105±25 
and 13.5±2.5, respectively) are characteristic properties 
of a starburst environment, where massive stars have had 
sufficient time to affect the isotopic composition of their 
surrounding ISM. The ^^C/^^C ratio is about twice as 
large as in the central region of the Milky Way, suggesting 
the presence of inflowing gas. 

Much work still remains to be done. Determining iso- 
tope ratios of sulfur in the nuclear environment of star- 
burst galaxies offers a unique chance to study the isotopic 
composition of sulfur ejecta from massive stars undergo- 
ing oxygen burning. Obtaining a more complete view onto 
the ratio in NGC 4945 and estimates of^^S/^^S in 

NGC 253, M82 and perhaps Arp220 would thus provide 
important insights into aspects of stellar nucleosynthesis. 
Detecting ^-^CN and obtaining a definite detection of ^'^CS 
would further constrain the ^^C/^^C and ^^0/^*0 ratios 
in the nuclear region of NGC 4945 and would thus also be 
fundamental for a better understanding of the chemical 
and physical conditions in a starburst environment. 

In spite of the large number of detections, our line sur- 
vey is not complete. While it appears that the chemical 
composition of the gas in the central part of NGC 4945 
is more resembling NGC 253 than M82, three important 
tracers have not yet been convincingly measured: thermal 
SiO, CH3CN and HCO. Furthermore, NO could be an 
interesting tracer for X-ray related chemistry. If our find- 
ing of similar chemical composition is valid for NGC 4945 
and NGC 253, SiO and CH3CN should be detectable in 
NGC 4945, while HCO should be weak. 

Our '^'^CO J=3-2 line profile is, with respect to inten- 
sity and lineshape, incompatible with the other CO tran- 
sitions observed. This also holds for the lineshape of the 
C3H2 2i2-loi transition with respect to the other C3H2 
line profiles. Another measurement with careful pointing 
is thus needed to obtain better agreement or to confirm 
the inconsistencies. 



By far the most important lack of observational data 
refers, however, to the absence of interferometric maps 
in NGC 4945. Molecular measurements are so far con- 
strained to linear scales of several 100 pc. It is still not 
known why some species show an almost Gaussian line 
profile (e.g. H2CO), while others show two pronounced 
peaks at ~450 and 710kms~^ (e.g. CN). High angular 
resolution maps could determine the location, kinematics 
and spatial distribution of the distinct velocity compo- 
nents. Where is the region with HNC/HCN>1 line inten- 
sity ratios located? Will high resolution HNC and HCN 
line maps provide important constraints for chemical mod- 
els of the starburst environment and the role of the AGN? 
And what is the distribution of the HNC/HN"C line in- 
tensity ratio that allows us to estimate optical depths? 
Can SiO, HCO and H^^CO+ be disentangled from each 
other as it was possible in the case of M 82? In the near fu- 
ture, interferometric maps of ATCA (Australian Telescope 
Compact Array) and SMA (SubMillimeter Array) will be 
available and single-dish telescopes like APEX (Atacama 
Pathfinder Experiment) will provide data on highly ex- 
cited spectral lines that are crucial to identify the major 
heating source(s). The goal of this paper is to provide a 
reliable observational basis for these measurements. 
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